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Synopsis 


This thesis introduces now nonstationary subdivision schemes for designing curves and 
surfaces. Any subdivision scheme for generation of curves consists of the following. 
Initially a set of points called control points such as 

po := : ae Z}, s = 1, 2 or 3 


is known. Any subdivision scheme {5a J iteratively generates a new set of control 
points 

:= € r : n € Z} 


at the itth level by a subdivision rule: 






The set : i € Z) of coefficients is called the mask at ^th level of the 

subdivision scheme. If tlu^ set a^*^) is such that its elements depend upon k, then the 
associated subdivision schouK^ is called a uoastationary subdivision scheme. 

Since the scluiine is applied comi)()n(mtwisc, one restricts the study to control points 
on E. Therefore, the initial data {p^ e R, a e Z} as well as the scheme {S.,} are 
considered on R only. 

The thesis consists of five chapters. The first chapter presents a brief review of sub- 
division schemes and some basic mathernatical preliminaries. In Chapters 2 to 4 three 
different schemes for designing curves are presented. For each scheme a subdivision 
algorithm .is constructi'd, its (•.onw'rgcncc', analysis is studied and some applications and 
advantages are outlined. Convergence of all tire schemes is shown using a theorem due 
to Dyn and Levin in [35]. In Chapter 5, a nonstationary subdivision scheme is pre- 
sented for generating smooth surfaces where the data points are not necessarily given 
on a rectangular grid. 



In Chapter 2, we present a subdivision scheme for the generalized Bernstein Bezier 
curves [38]. The scheme is a generalization of Lane and Riesenfeld algorithm for clas- 
sical Bezier curves [49]. As particular cases one can obtain subdivision schemes for 
classical as well as trigonometric Bezier curves [48]. We also point out the advantages 
of using the subdivision scheme for designing curves as well as surfaces. The algorithm 
is given as follows. 

Algorithm 0.1. Starting from q^, i = 0,1, ... ,n the set of control points qf, i = 
0,1, ... ,2n are obtained as follows: 

Letpl = (?f, i = 0,l,...,n. 

For k = 1, . . . ,n, set 

Pf = boV ) Pti' + (-.(0 Pb'. i = (; + 1. . . . n. 

Define 

1 j Pit "i 0, 1, . . . , n, 

4 o • 

[p;" 2 = n,n + l,...,2n. 

In the a])ove algorithm /.)()(/) = and bi{l) = where d(x) is the solution of 
an initial value problem [38]. Besides the algorithm we also give some useful properties 
of the function d{x) which generates the space of generalized Bernstein polynomials 
and some results conccnaing geiierHlized convex hull, called 2?-convex hull associated 
with d{x). 

There is extensive work involving trigonometric spline as a tool for designing curves. 
Successful efforts have been made in [48] to establish algorithms for trigonometric 
splines analogous to algorithms for ])olynoniial splines of one variable. There is an algo- 
rithm in literature for trigonometric splines [48] analogous to Oslo algorithm. However, 
there is no algorithm' for trigonometric splines corresj)onding to Lane- Riesenfeld algo- 
rithm of polynomial spline. In Chapter 3 we introduce an analog of Lane-Riesenfeld 
algorithm for the trigonometric spline with uniform knots. It turns out to be a non- 
stationary binary subdivision scheme. 

Algorithm 0.2. Lei n > 1 and n — ‘2k — 1 or n = 2k for some k > 1. Starting from 
qf, i — 0, 1, . . . , 7i the first level conlml points q\, i — n — 1, . . . , 2rn + 1 are obtained as 


follows: Set 



qI 

flil (nP 

a{l) Wi- 


+ 9°). 


1 = 0,l,...,m if n = 2k — 1 

2 = 1, 2, . . . , m if n = 2k 


Pii+l Qi > 


i = 0, 1, . . . , m (n = 2k orn — 2k — 1). 


For j = 2, 3, . . . /c, define 

w{2j - 3, 1) {w{2j - 2, 1) piZ 2 + Pill + w(2j - 2, 1) 

j for n = 2A; — 1, and i = 2j — 2, 2j — 1, . . . , 2m + 1. 

w{2j - 2; 1) {w{2j - 1; 1) jfizi + + '^(2i - 1; 0 

for n = 2k, and i = 2j — 1, 2j, , 2m + 1. 

Set Qi = Pi, z = n — 1, . . . , 2m + 1. 

In the above algorithm w{j]l) = • The basic mathematical contribution 

consists of an introduction of a convolution formula for trigonometric B-splines. The 
above scheme heli)s iti the exact reconstruction of cosine, sine functions and circles. 

In Chapter 4, we introduce a nonstationary subdivision scheme whose limit curve 
interpolates a given set of data points. 

Algorithm 0.3. Given the control points {pP € K, z = —2, —1, . . . , rz + 2} the control 
points , i = —2, — 1 , 0 ,..., 2^‘ ' + 1 } at level k + 1 are obtained by the folloxving 

recursive relation: 


Ihi 


it 


■1 < i < 2^n + 1 


/4i+l — -"^^’kltl T (2 '<‘>k)pi + (2 + '^‘’k)Pi+l 




■1 < z < 2'=n 


where Wk 




2sin(J^)sin(p^-)' 


This scheme is a generalization of a well known four point subdivision scheme of 
Dyn and Levin [32]. We establish that our scheme is asymptotically equivalent to the 
four point schenu! of Dyu and Levin [35]. One of the important features of this scheme 
is that it niproduces functions 1, cos(aa:), sin(o;.T) and their linear combinations. In 
particular if we clioos(! vertices of a regular zt-gon as initial control points, then the 
limit curve turns out to be the original circle passing through the vertices of the regular 



n-gon. We also show that this scheme can be used to approximate srnootli functions 
quadratically. 

The above schemes for curves arc generalized in a natural way to nonstationary 
subdivision schemes for tensor product surfaces. These schemes require that the data 
set should lie in a rectangular grid. However, quite often the data set does not lie on 
a rectangular grid. There are several stationary subdivision schemes like ones given 
by Doo-Sabin [26] and Catrnull-Clark [16] which generate surfaces from such data. 
However, there are only few results available dealing with nonstationary subdivision 
schemes generating surfaces for such data except the work of Dyn and Levin [35] which 
deals with tensor product exponential splines and exponential box splines. Moreover, 
formulating nonstationary schemes for arbitrary data and their convergence analysis 
is a difficult task since the imusks keep changing in each iteration. In Chapter 5, we 
present a nonstationary subdivision scheme generating surfaces for data prescribed on 
an arbitrary topology. The scheme can be viewed as a generalization of the subdivision 
scheme for the trigonometric splines introduced in Chapter 3 to cases with arbitrary 
surface data. The schenu; is formulat(‘d along the line of the Doo-Sabin and Catmull- 
Clark schemes and convergence analysis follows the mechanism of U. Reif [63] and Jorg 
Peters [59]. 


viii 
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Chapter 1 


Introduction 


1.1 A Brief History 

Computer Aided Geometric Design or CAGD in short evolved during 1960s as a result 
of development of digital computers and simultaneous development of software for 
designing curves and surfaces in designing industry. These new tools to a large extent 
replaced expensive methods of building models and templates. Therefore the focus 
shifted to find computer compatible forms of curves and surfaces by replacing the 
master model by a mathematical model, a computerized form of 3D object. This 
resulted in defining a surface as a set of triplets in and a smooth injective map from 
R^toRT 

The parametric curves and surfaces which were till then exclusively studied by 
differential geometry [13] were subjected to numerical studies on the basis of their 
numerical construction ushering a new area of study called Computer Aided Geometric 
Design formally so designated in 1974 [5]. 

In surface reconstruction, one is provided with a large number of point coordinates. 
This information is obtained from the drawing board or from the physical model. There 
exist several methods to generate surfaces from this set of discrete points. Some of the 
important methods like the Bezier methods, B-spline methods, and subdivision meth- 
ods are briefly outlined below. 

Bezier methods: 

The Bezier methods of generating curves and surfaces were developed independently by 
P. de Casteljau about 1959 at Citroen [14, 15] and by P. Bezier about 1967 at Renault 
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[6, 7, 8, 9]. The work of Bezier which was titled UNISURF was published earlier than 
the Citroen technical report and hence the theory now bears the name of Bezier. 

Bezier curves and surfaces j)rovide a powerful tool in curves and surfaces design: 
They are affinely invariant, stable and have linear precision property. Moreover the 
Bezier curve has variation diminishing property. 

In theory of Bezier curves barycentric coordinates and Bernstein polynomials play 
a key role. This encouraged others to define barycentric like coordinates, and define 
Bernstein polynomials and Bezier curves on these newly defined objects. Recently, 
Alfeld, Nearntu and Sdiumaker [1, 2] introduced barycentric like coordinates for cir- 
cular and s[)}ierical ooonliuates. This helped to introduce circular Bernstein-Bezier or 
CBB c:nrves and si>herica! Bi'rustein-Bczier or SBB patches. The CBB curve is also 
called as trigonometric Bernstein-Bezier or TBB curve by Koch et.al.[48]. Before them 
Brown and Worsey [12] had tried to define barycentric coordinates on sphere. Recently 
Gonsor and Nearntu [38] presented a unified Bernstein-Bezier theory via the concept 
of a polar form. They consider a class of functions (here we call such functions gen- 
eralized Bernstein-Bezier or GBB polynomials) which include the classical Bernstein 
Bezier polynomials, trigonometric Bernstein Bezier polynomials as special cases. The 
associated barycentric: coordinates are defined as solution of an initial value problem 
of a second order differential equation. The concept of Bernstein-Bezier polynomial is 
now very broad which contains polynomials as well as non polynomials. 

B-spline methods: 

Study of splines started with I. J. Schoenberg in 1946 as univariate piecewise polynomial 
functions [67, 68] and has been extensively studied in the theory of approximations since 
then. 

Apparently the earliest use of spline in CAGD was about 1963 by J. Ferguson [37] 
in Boeing Company. He used cubic splines to interpolate a given data. At about 
same time C. de‘ Boor and W. Gordon studied those curves at General Motors. Use 
of splines became popular in CAGD after the discovery of compactly supported basis 
function called B-splines. The name B-spline was first used by Schoenberg in 1967 
[70]. Schoenberg first dealt with B-spline on equispaced knots and defined them via a 
convolution formula [67]. This convolution formula helped to find an efficient algorithm, 
known as Ijuiu) and Ricsenlcld algoritluii for univarite Hplinci' curve. Gordon [39] and 
Riesenfeld [65] first used the B-splines in the context of parametric curve design and 
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showed that they are proper generalization of Bezier curves. A B-spline can also be 
expressed as a piecewise Bezier form [42]. 

The spline curves have most of the properties of Bezier curves; They are affinely 
invariant, variation diminishing and possess linear precision property. In addition to 
that B-spline curves have local control property. 

Since the discovery of splines many successful attempts have been made to gener- 
alize B-splines, involving class of functions which are not polynomials such as rational 
splines, trigonometric splines, Chebychev splines and L-splines etc. A comprehensive 
theory on the.se uon-i)(.'lynomial spliiuis can be found in the book [71] by Schurnaker. 
The contemporary meaning of spline is very broad and it includes general piecewise 
functions rather than piecewise polynomials only. 

Recently trigonometric splines were given special attention for designing curves. 
Ti'igonometric splines wore introduced by I. J. Schoenberg in [G9], and have been exten- 
sively studied in the literature [47, 48, 54, 55, 71, 72, 74]. It has been observed that 
most features of classical polynomial splines carry over to trigonometric splines, in par- 
ticular, for example, it is shown that any trigonometric spline is a linear combination 
of trigonometric B-splines [55, 71]. Besides, a recurrence relation for trigonometric 
B-spline has been established in [55]. In [48], the notion of control curve for a trigono- 
metric spline has been introduced and studied, besides the convex hull property and 
variation diminishing property of trigonometric splines arc also established in the spirit 
of the results available for polynomial splines. / 

The simplest procedure to obtain two dimensional splines is to use tensor prod- 
uct for rectangular regions. The tensor product splines inherits all the properties of 
univariate splines except variation diminishing property. A more general approach 
to the construction of multivariate splines consists of simplex splines [10], polyhedral 
splines [22] and box splines [11]. Multivariate splines are extensively studied in [20, 57]. 

Subdivision methods; 

Subdivision methods consist of classes of numerically stable easily implementable al- 
gorithms for the generation of parametrized curves and surfaces. These algorithms are 
sometimes referred to as corner cutting algorithms. Essentially, a subdivision method 
describes a continous shape as a limit of a sequence of discrete shapes described by 
so called control points. Tluise discrete shapes are obtained by linear transformations 
based on splittitig and aveuaging. 
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Subdivision methods are also faster and easier methods to construct Bezier curves 
and surfaces as well as spline curves and surfaces. Moreover subdivision schemes are 
also used to generate surfaces from a data set associated with an arbitrary topology 
in a very elhcient way unlike other available methods. They are also intimately linked 
to multiresolution structures and are used widely in other areas of computer graph- 
ics such as constructing wavelets and fractals, animation and rendering etc. Recently 
subdivision schemes are applied to solve problems in variational calculus [45, 75], ordi- 
nary differential equations [70, 78] and to approximate functions and convex compact 
sets [30]. 

The general set up for design of surfaces via the subdivision method is the fol- 
lowing: The designer starts from a finite set of points, called control points and a 
connectivity rule between them. The points are arranged according to the univari- 
ate index set Pi, P2, ...,Pn € (or E^). Then more and more control points are 
generated by a subdivision rule. At each iteration level a control polygon is iden- 
tified with the control points. If the subdivision rule is suitably chosen, the sequence 
of control polygons converge to a smooth limit function in the sense that the limit 
function has the similar shape as the control polygon and is very close to them after 
certain iterations. The curve can also be made to pass through the original control 
points. In this case, the schemes are called interpolatory. 

The general set up for surface designing in the subdivision method is the following: 
We start with a finite set of control points in and a connectivity rule among 
them. The connectivity rule associates the control points to vertices of edges, faces 
and hence associates them to a control polyhedron. The oder of a face is the 
number of vertices constituting the face. The order of a vertex is the number of 
faces associated with the vertex. The control polyhedron is said to have a regular 
topology if the control points have an indexing 

p,^ G R‘ , a 6 

and order of all the faces, order of all the interior vertices and order of all the boundary 
vertices are same respectively. The regular topology contains a tensor product topology 
as well as a triangular t()|)()logy.' Tlien a subdivision rule is used iteratively for 
generating new sets of control points and corresponding control polyhedrons (see Fig 
1.1). The rule is such that the seriuence of control polyhedrons converges to a smooth 
limit surface. 
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Figure 1.1: Catmull and Clark subdivision scheme 

In a subdivision scheme, if the subdivision rule used is same at all levels of iterations 
hen the subdivision scheme is called a stationary subdivision scheme otherwise it is 
called a nonstationary subdivision scheme. 

There are various stationary subdivision schemes in the literature. One of the 
earliest stationary subdivision schemes for curves was studied by G. de Rham [64]. 
In 1974 G. Chaikin presented a subdivision scheme for quadratic spline curve [19]. 
Lane and Riesenfeld generali'/ed tlie Chaikin algorithm and presented a subdivision 
scheme for polynomial spline cur ye of arbitrary degree [49]. The latter subdivision 
scheme is considered as one of the important stationary subdivision schemes in CAGD. 
Another important stationary subdivision scheme for design of curves is the 4-point 
interpolatory scheme by N.Dyn [32] (see also [27]). ' 

Stationary subdivision schemes for curves are extended to generate tensor product 
surfaces in a natural way if the control points are arranged in a rectangular topology. 
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For triangular topology there are also subdivision schemes such as line averaging algo- 
rithms [24], subdivision schemes for box splines [23, 43] and butterfly scheme [33]. All 
the above subdivision schemes are studied by a unified theory as given in the work of 
Cavaretta and Micchelli [17] and Dyn [29](see also [18, 31, 56]). 

Stationary subdivision schemes for arbitrary topologies were first studied in 1978 
by Catmull and Clark [IG], Doo and Sabin [2G], and later by Loop [52]. Ball and Storry 
[3, 4] analyzed these schemes using a matrix approach. In 1995, U. Reif presented a 
theory for these schemes [62, 63]. He provided a necessary and sufficient condition for 
the convergence of such schemes ensuring tangent plane continuity of the limit surface. 
His study was based on the nature of eigenvalues of the associated local subdivision 
matrix. After this, a stream of new schemes like mid-edge subdivision scheme [58], 
Kobbelt subdivision scheme [44, 4C, 73] , tagged subdivision scheme [41] and combined 
subdivision scheme [50, 51] were introduced. At the same time new approaches to 
analyze these schemes were also introduced [59, 60, 61]. 

Recently several nonstationary subdivision schemes for curves have been introduced. 
The first nonstationary subdivision scheme for curves was proposed by Rvachev [66]. 
In [35], Dyn and Levin presented a general theory for a special class of nonstationary 
subdivision schemes, called asymptotically equivalent subdivision schemes, where, as an 
example they presented nonstationary subdivision schemes for exponential splines and 
exponential box splines. A nonstationary scheme generating circles is also presented 
by Dyn in [34]. Nonstationary subdivision schemes are also recently used by .1. Warren 
to solve problems in variational calculus [75] and by H. Weimer to solve ordinary 
differential equations [77, 78] and physical problems in fluid dynamics [76]. 

1.2 Contribution of the Thesis 

In Chapter 2 of the thesis, we ■ present a subdivision scheme for the generalized 
Bernstein-Bezier curves. This scheme is a generalization of Lane and Riesenfeld al- 
gorithm for classical Bezier curves [49], As particular cases of the above scheme, one 
gets the subdivision schemes for classical as well as trigonometric Bezier curves. We 
also highlight certain- advantages of using this subdivision scheme for designing curves 
as well as surfaces. Bcisides tin; algorithiu w(! also (hiscribe some useful pr()i)erties of 
the function d{x) which generates the space of generalized Bernstein-Bezier polynomi- 
als. Besides, some results dealing with generalized convex luill, called, X>-convex hull 
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associated with d{x) is also discussed. 

Recently researchers working in the area of designing of curves and surfaces have 
started looking trigonometric spline as a tool for designing purposes. Successful ef- 
forts have been made [48] to establish algorithms for trigonometric splines analogous 
to algorithms available for polynomial spline functions of a single variable. There is 
an algorithm in literature for trigonometric splines [48] generalizing Oslo algorithm. 
However, there is no subdivision algorithm for trigonometric splines corresponding to 
Lane-Riesenfeld algorithm of polynomial splines. In Chapter 3 we introduce an analog 
of Lane-Riesenfeld algorithm for the trigonometric splines with uniform knots. A non- 
stationary binary subdivision scheme is introduced to generalize this algorithm. The 
basic methodological contribution in this Chapter is the introduction of a convolution 
formula for trigonometric splines and showing the probability of exact reconstruction 
of parametric curves (cos (x), sin (x)) or circles. 

In Chapter 4, we introduce a nonstationary subdivision scheme whose limit curve 
interpolates a given set of data points. This scheme is a generalization of a well known 
four point subdivision scheme of Dyn and Levin [32]. We establish that our scheme is 
asymptotically equivalent to the four point scheme of Dyn and Levin [35]. One of the 
important features of this scheme is that it reproduces functions 1, cos(x), sin(x) and 
their linear combinations. In particular if we choose vertices of a regular n-gon as initial 
control points, then the limit curve is the original circle passing through the vertices of 
the regular n-gon. We also show that tins scheme can be used to approximate smooth 
functions effectively. 

There is not much literature available dealing with nonstationary subdivision schemes 
for the generation of surfaces except the work of Dyn and Levin [35] which presents 
nonstationary schemes concerning tensor product exponential splines and exponential 
box splines. In Chapter 5 of the thesis, scheme for trigonometric splines introduced 
in Chapter 3 is generalized into a nonstationary subdivision scheme generating surfaces 
from tensor product topologies as well as arbitrary topologies. The subdivision scheme 
is in some sense a generalization of the Doo-Sabin scheme and the Catrnull-Clark 
scheme. The limit surface retains same shape as the limit surfaces of these stationary 
schemes, but the sizes of the limit surfaces can be varied by changing certain design 
parameter or functions. 
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1.3 Mathematical Preliminaries 

In this section we briefly outline some relevant definitions and results on subdivision 
schemes. First we consider subdivision schemes concerning a regular topology. We 
mainly concentrate our discussion on nonstationary subdivision schemes. 

Given a set of control points 

:= {Pa € : a G Z^}, d = 1 or 2, s = 1, 2 or 3 

at level 0, a subdivision scheme generates a new set of control points 

:= (pf. e r : fv e Z''} 

at the kth. level by a subdivision rule: 

pS = • ■ • S,P°)s = Y, 4-2= rt', /3 S Z". (1.1) 

The set := : i € Z‘^} of coefficients is called the mask at kth level of the 

subdivision scheme. If the mask is independent of k then the scheme is called stationary 
otherwise it is called uonstationary. To each subdivision scheme {S,,*.} defined by the 
masks we a.ssign kth level characteristic multinomial (polynomial) 

^ A: > 1 

called the A;th level characteristic polynomial of the algorithm. 

As the schemes of the type {-Sa*.} are applied componentwise, it is enough to study 
subdivision schemes for the initial data P” = {p° : a 6 Z*^}. In these subdivision 
schemes we associate the points k 6 Z'^ to the location values a G Z*^. 

Definition 1.1. [35] A subdivision scheme {Pa,. } is said to be if for every initial 
data P° G l°° there exists a limit function / G such that 

lim sup Ip* - /(2-*o)l = 0 (1.2) 

and / ^ 0 for some initial data P°. 

Definition 1.2. [35] A subdivision scheme {Pojt} is called stable if there exists a 
constant M such that 


Vm.;c€N. 


(1.3) 
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The limit function F obtained by the application of the scheme on the initial 

data 

[O, 

is called the basic limit function of Here 0 refers to the zero element of Z‘^. 

The following Theorem, which is immediate from Theorem 13 and Lemma 14 of [35] 
is one of the basic convergence characterizations for general nonstationary subdivision 
schemes. 

Theorem 1.1. Let a sequence of continuous Junctions {ipm ■ > 0 satisfy 

(a) lim j) = l uniformly in W'' 

m-voo ^ 
ieZ‘' 

(b) ||'0m||oo < B for some constant B 

(c) 'ijjrn, m >0 are compactly supported. 

If the non- stationary scheme is and L is its limit function for the initial data 

P° then the seq'uence of functions {f^ : — >■ R} given by 

f^{x) = ^p^M2^x-i), . (1.5) 

converges uniformly to L on the support of L. 

Conversely, if a sequence of functions {ipm € C'(R‘^)} satisfies the stability condition: 

''ll»ll» S -‘)llco- 9 e '”(2"). s (l-O) 

ie'£‘‘ 

for some K E and the sequence of functions {/^} converges uniformly to a contin- 
uous function /°° on compact subsets of R^, then the BSS {Paj,} is a convergent 
subdivision scheme. Moreover {Pa*,} is C’^ if f°° € C''”(R). 

In particular if d = 1 and : R R, m > 0 are given by 

fpmit) = (1 - \t\)+, (1.7) 

whore ( denotes the i)lus fuii(;tion or the truncated error function, tlien xjjm is com- 
pactly ; |) orted and clearly satisfies (a), (6), (c) of (1.4). The functions {/*} represents' 
the S 3 ( | fij|e of control polygons joining the points of P^. An obvious consequence of 
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the above Theorem is that if is (7° then the sequence of control polygons {/*} 

converges uniformly to L. In practice, taking infinite number of iterations is not possi- 
ble, therefore, the import of the Theorem is that the polygon after certain number 
(say k) of iterations can be taken as an approximate representation of L. In case d= 2, 
analogous polygons are used for representing the limit surface [40]. 

Moreover, »>■ > 0 satisfy (l.G). Therefore if the sequence {/*} associated with 
this converges uniformly then the subdivision scheme {S'aj.} converges. This is one 
way of showing convergence of a subdivision scheme. 

Let Ed denote the set of extreme points of a unit d-cube. In particular Ei = {0, 1} 
and E '2 = {(0,0), (0, 1), (1,0), (1, ])}. Let us define the matrix norm corresponding to 
the scheme {Sa^ } by 

L == - (1-8) 

0ez<‘ 

and define the concept of asymptotic equivalence as follows: 

Definition 1.3. [35] Two subdivision schemes {^a*,} and {Sb^} are asymptotically 
equivalent if 

00 

\\^0.k ~ ^bic lloo 

k=l 

Moreover if {Sb^} is a stationary scheme then {5a*,} is called an asymptotically sta- 
tionary scheme. 

Let us define the support of a mask as the set {a G ^ 0}. As the 

following Theorems show the convergence of two asymptotically equivalent schemes 
are closely related. 

Theorem 1.2. [35] 

(a) Let {5a^} and {56*^} be two asymptotically equivalent schemes. Then {5a;(,} is 
and stable if and only if {5(,*,} is and stable. 

(b) //{5a*,} is asymptotically equivalent to a stationary scheme {5a} with a finite 
mask then {5a*,} is and stable. 

Theorem 1.3. [35] Let {5a*,} and {5„} be two a.symptotically equivalent subdivision 
schemes having finite m.asks of same .support. If {Sa] is and 

(XJ 

k=0 


( 1 . 10 ) 
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then the. non-stntionary BSS {‘S',,;,. } *.s' C"'. 

The following examples provide some iionstatioriary subdivision schemes for curves. 

Examples 1.1. (Exponential B-spline) [35] The non-stationary subdivision scheme 
{Sa^} having kt\i level characteristic polynomial 

m 

i-1 

is a convergent scheme. Moreover the basic limit function F is the exponential 

B-spline with integer knots, spanned by j = I, . . . ,m} and supported on (0,m). 

Examples 1.2. (Circle generating scheme) [34] The scheme is represented by the 
characteristic polynomial 

with ak = 2~^~^aQ, 0 < ao < tc luis the limit function in ^'(R). Moreover, if the 
initial control polygon is a regular n-gon and we choose cxq — then the limit curve 
is the circle inscribed in that regular n-gon. 

The following subdivision scheme is an interpolatory subdivision scheme for design 
of curves. 

Examples 1.3. (Four Point Scheme) [27, 32] The scheme is defined by the mask 

a := {o _3 = 03 = —w, a_i = Oi = ^ -f w, Oq = 1}, Oj = 0, f —3, 3, 0, —1, 1. 

It has been known that the scheme converges for |ro| < 1/2 and the limit function 
belongs to C'^(R) for the range 0 < w < 0.154 [28, 29]. Moreover the associated 
characteristic polynomial a{z) can be written in factorized form as given below [29]: 

a{z) = z~^{l + z)'^(^z'^ — 'w{z — 1)^(1 + z"^)). 

jL 

The subdivision schemes for curves are easily generalized into schemes for surfaces 
with regular topology. In case, data set is prescribed on grids with irregular topologies 
uniform sul)divisioa rule cannot be applied to all the faces and control points. Therefore 
it is of grtsat importance to develop subdivision schemes applicable to such situations. 
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Doo-Sabin [26] and Catmull-Clark [16] (see Figure 1.1) first developed such schemes. 
We describe below the general procedure for these two schemes. 

It is enough to consider the set of initial control points as vertices of a polyhedron 
which has a face of order n sorrounded by a layer of regular meshes. Let denote 
the set of control points at A:th level. Then the set P^'+’ is obtained from the set 
by a matrix operation [63]; 

= MP^. (1.12) 

The matrix M is called the subdivision matrix associated with the scheme. In the 
analysis some results from the theory of circulant matrices are used. 

Definition 1.4. (Circulant Matrix) [25] The block matrix 

^ j4q A\ • ■ ■ 

An-\ A{) ••• /1„_2 

A := 

Y y4l A 2 • ■ • Aq 

where Ai, i ~ 0, 1, . . . ,n - 1 cire matrices of order {m,m), is called a block circulant 
matrix of order (n, ni) and is denoted l)y A = bciTa{Ai), A\, - • ■ , An-i)- 

Denoting by u = := exp{i^). Define the Fourier matrix Fn as 

1 1 1 

1 u '(P 

1 

where F* is the conjugate transpose of a matrix P. It is known [25] that 

(z)P„ = P„L (n)p-i = P„*, (m)(P:)-i=P„. 

For A = bcirc(ylo, vli, . . . , A„_i), the m x m piatrices A*;, /c = 0, 1, . . . , n - 1 are 
defined by 

n-l 

j=Q 





The following result holds [25] 
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Theorem 1.4. (Diagonalization of a block circulant matrix) A block circulant 
matrix A = bcirc{Ao, , An-i) has the diagonalization 

^ = (F„ X FmTdiagiFraAoF:;,, • ■ • , x F^). 

The following is a corollary of the above diagonalization result. 

Corollary 1.5. [25] Let A be an eigenvalue of a particular matrix Aj, then it is also 
an eigenvalue of An- j as well. Moreover, ifv is an eigenvector of Aj corresponding to 
A then the vector v defined by 


V := [f;, •, (1-14) 

and V, the complex conjugate of v are also eigenvectors of A with eigenvalue X, where 
T denotes the block tramposc. 


It is known that in case of Doo-Sabin and Catmull-Clark schemes the matrix M of 
(1.12) is a block circulant matrix [63]. This class of schemes have been recently studied 
in a unified way by U. Reif [63]. He has shown that the Doo-Sabin class of subdivision 
schemes converges and the limit surface has a continuous tangent plane everywhere. 

Briefly, one variant of Doo-Sabin scheme is described by the following. Let F be 
a face of with n vertices A\,. . . , An- Then the scheme defines the new vertex Oj 
corresponding to AljS by the subdivision rule: 


4n + 2 7). -|-2, ^ A A \ 

m = — r Ai H + -^i-l) + '^{Ai+2 + • • • + j4i+n-l)- 

an 


8n 


8n 


A comprchensivo collection of all available subdivision schemes for arbitrary topol- 
ogy and their analysis is contained in [73]. 




Chapter 2 

Generalized Bezier curves 


In this chapter a biliary subdivision scheme for evaluation of the generalized Bernstein 
Bezier curves is presented. This scheme is a generalization of the Lane-Riesenfeld 
algorithm for polynomial Bezier curve. Moreover in this chapter we present a procedure 
to find the P-coiivex Imll of a sub.set of E^. 

2.1 Generalized Bernstein-Bezier Representation 

In this section we recall some basic concepts of the unified Bernstein-Bezier theory [38]. 

Let P be the two dimensional null space of a second order constant coefficient 
differential operator L of the form 

L-.= 7,(lieR. (2.1) 

Let d E 'D be the unique solution of the initial value problem 

Ld = 0, d(()) = 0, Dd{{)) = 1. (2.2) 

It follows that V = span{d(- — t) : t e E}. Since V is two dimensional, for any two a 
and b such that d{a — b) ^ 0 (and therefore d{b — a) ^ 0) we have 

V = span{d(- — a), d{- — /;)}. 

For n > 0, let P„ ;= span{d"(- - t) : t € R}. The elements of T>n are called V- 
polynomials of degree n. 

The following result is proved in [38] . 
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, , _ , ,, thni n <h and d(b - a) ^ 0. Then the functions 

Theorem 2.1. lei a, 5 € R be such that a k j t- 


n 


Bnx-,a,b) := , 


d(^J 


a: e R, z = 0, 1, . . . , n (2.3) 


%= I 


form a basis for Vn . 

, . . . n 1 n the Bernstein basis polynomials of degree 

We call the functions Sf , z = 0, 1, • ■ • . 

n coressponding to V or simply ^-polynomials. 

^ , r ^ ,!,>,> T) = rin, where n.„ is the space of algebraic 

In case 7 = 0 ^ = 0, ‘ n- 

^ ^ , i-hor hand if 7 = 0 and <5 = 1 then we have Vn = 7n, 

polynomials of degree n. On the other liana j 

where 

fsp,an{l, «n{2x),»s{2*)..«(4x); cos{4x ).. . . , sin{nx), co.(nx)}, 

n ev(!n, (2.4) 

spa,.{xm{.x), .™(x), x«(3x), co»(3x), .... sin{nx), cos{nx)], 

n odd, 

r * • rnlvnonmls of order n + 1 [1, 48]. Clearly when 

is the usual space of tncjonomeino poiynu, 

„ , „r, • _ n 1 n are the classical Bernstein basis poly- 

Vn = n„ the ^-polynomials 5,", z - U, -L. ■ • • ' _ ^ 

nomials [30] and in caac C„ = %. iba B-polynomials are the circular Bernstern basrs 
polynomials ID. or trigo.romctric Bernstein basis polynomials |48|. 

T . rr, ,x.W that d(b - a) ^ 0. Then for any Pi G R, 

Definition 2.1. Let a,b € R such tnar, 

z = 0,1 ,..., n we call 


2=0 


(2.5) 


a generalized Bern.stein-Bfaier or GBB-polynomi,al of degree n. 
The points {(Zi> c.) : t = 0, 1, ■ - 'tt}. where 

b~a 


^i = a + ^- 


n 


i = 0, l,...,n, 


(2.6) 


are called control points of / . 

Let us define £„ := span{ci(n ■ -t) - te R}- 

Definition 2.2. (Control curve): , , . ^ . 

. /o r'l Let r/i, * = 1, 2, ..., n, be the unique functions 
Let /(.r), X e [a, b] be defined as m (2.oj- 
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from which interpolate the control points i = 1, 2 , . . . , n. More 

precisely, define 


9i{t) ■■= 


d{n{t - ^i)) d{n{t - 

rf(nte-, - {i)) d(n(ii - 5i_0) 


(2.7) 


We call the curve consisting of the i)ieces Gi := {{t,gi{t)), t € i = 1, . . . ,n, 

the control curve of /. 


Definition 2.3. (X>-convex set): 

A subset B of is called X>-convex of degree n, n > 2 if for any two points (ti,ci) 
and (^2, C2) of B such that d{t2 — ti) 0, the curve 


r f. d{n{x - t2)) d{n{x-ti)) . 
d{n{U-t2)) "'■^d(n(t2-ti)) 


t € [ti,t2]} 


which is the P-convex combination of degree n of (ii,ci) and (12,02) lies in B. 


Definition 2.4. (X>-convex hull): 

The D-convex bull of degree n of a subset B of R^ is the smallest X>-convex set of 
degree n containing B. 


Trigonometric convex hull of degree n which is introduced in [1, 48] is a special case 
of P-convex hull. 

If h{t) is a GBB-polynornial on [«,/)] with control points {(ti.Ci), i = 0,1 ,..., n} 
then it can be shown [38] that {(t, h(t)) : t E [a, 6])} lies in the P-coiivex hull of degree 
n of its control points {(tj, Cj) : i = 0, 1, . . . , n}. 


2.2 Properties of d(x) 

We now give some properties of d(x) which will be useful in the subsequent sections. 
We recall here that d(x) is tlie solution of the initial value problem (2.2). Let 
+ r7 + 5 = 0 be the characteristic equation of the differential equation Lf = 0, 
where L is the differential operator defined by (2.1). Let r\ and r2 be the roots of the 
characteristic eciuation. Depending upon the values of rq and r2 we have the following 
cases: 

(A) d(x) = = 2 - ^ if and r2 are real and distinct. 

(B) d(x) = xeC^^ = X if 7-^ = — -y. 
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(C) d{x) = e s\xi{j5x)/ jd if = ct + i/5 and r 2 — ct — ip. 

Note that d{x) > 0 for a: G (0, z) where z = oo for the cases (A) and (B) and z = Tij jd 
for the case (C). Moreover d{x) can have at most one zero in any interval I of length 
|/| < z. Therefore the function d{- — t), i 6 M can have at most one zero in any interval 
I of length |/| < z. 

Observe that d{nx) is the solution of the initial value problem: 

Lnf = 0; /(O) = 0 ; f (0) = n (2.8) 

where := + nyD + S and £„ is the null space of the operator Moreover 

d{n •) and its translations d(n • —t), t G E can have at most one zero in any interval 1 
of length |7| < A This leads to the following lemmas which we use in the sequel. 

Lemma 2 . 2 , Let n(x) G Then n(x) can have at most one root in any interval I 
of length |/| < z/n. 

Proof: Let a:o G R be such that u{xq) = 0 and u'(xo) = M for some M G R. Then 
u{x) is the unique solution of the initial value problem: 

Lnif) = 0; f{xo) = 0; f{xo) = M. 

Hence u{x) = ^d{n{x — a:o)). Therefore u{ic) 7 ^ 0 for \n{x — a;o)| < z: or |a: — a:o| < 
Thus the lemma follows. □ 



Figure 2.1: Two distinct solutions of £„ intersect each other only Once. 

Remark 2.3. Let a,h G R such that 0 < b — a < zjn. Let /i, /2 G £n such that 
/i 7 ^ /z- Let /i(xo) -- ibr some Xq G Then by Lemma 2.2 and the fact 

that /((xo) 7 ^ f^i^o), we have either (see Figure 2.1) 

f\{x) < / 2 (x), X G [a, xo) and fi(x) > / 2 (x), x G (xq, b] 



Chapter 2 


18 


or 


h{x) > f2{x), X e [a, xo) and fi{x) < f2{x), x e (a;o, 6]. 
Lemma 2.4. Let a,b eM. and 0 < b — a < Then the function 

_ sinWx - g)) smWb - x)) ^ 

s\n{p{b — a)) sin(/3(6 — a)) ’ 

has the rnaxinium at Hence inax{i7(x) : x E [n, b]} = 

Proof: We first note that ' 

'cos{P{x — a)) cos{P{b — x)) 


g'{x) = ft 


, X E [a, 6 ] 


H\\\{ft{b — a)) H\n{ft{b — a)) 

is a solution of the differential operator L := D"^ + ft"^. Therefore g' can have at most 
one zero in [a, 6]. Since g'{^) = 0, ^ is the only zero of g' in [a, b]. Noting that 
g(a) = g(b) = 1 and — is a strictly decreasing function in (0, |), we have 

, . X - a b - X , ^ , ,, 

Oix) > 7 + 7 = 1 Vo: € (a, b). 

b — a b — a 

Therefore, g has maximum equal to ^ at x = □ 


In an analogous manner, one proves using strictly increasing property of in 

(0, oo), the following Icmima. 

Lemma 2.5. Let a,b E and i) < b — a < oo. Then the function 

sinh(^(x — o)) sinh(/3(6 — x)) 


g{x) 


+ 


sinh(,0(6 — a)) sinh(/3(f) — a)) 


-, X E [a, 6] 


has the minimum at Hence min{g(x) : x 6 [a, 6]} = 


cosh(^^)- 


In the following proposition we give some properties of d{x) which will be used in 
finding certain estimates in the next sections. 

Proposition 2.6. For 0 < I < z, we have: 

>0 V,;6(0,il. 

(ii) Either M V.T 6 (0, (] or V.x e (0. i], 

(ni) Bilker > 1 Vl € (0, (] or < 1 Vx € (t), i|. 
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(iv) There exists a constant Ci depending only upon 7, S and I such that for all 
X € [— /, /] we have |fi(x)| < Cila:!. 

Proof: The values of for the cases (A), (B) and (C) are and 

respoctivoly. Similarly the values of for the case (A), (B) and (C) 
are ie-^1^/2 and 6““^/^ respectively. This proves (i) and (ii). 

The values of for the cases (A), (B) and (C) are 

cosh(7a;/4) and respectively. This proves {in). 

The proof of (iv) follows from the mean value theorem for the function d{x). □ 


The following lemma is also needed in the sequel. 

Lemma 2 . 7 . Let 0 < I < z and 0 < v — u < 1 . Then there exists a constant C4 
depending only upon 7, 5 and I such that for all x G [u, t;] 


h{x) :■ 


Moreover, in case 7 = 0 


1 

H 

1 

1 

1 

1 

1 

1 

p 

1 

1 


1 

1 

1 


< C4 (r; — w) . 


h{x) < C4{v — uY, X G [u,u] 

g different 
and h2{x) 

h\{x) < _ i| and h2{x) < — 1 


Proof: We prove the lemma by considering different cases ((B), (C) and (A) of this 
section) of d{x). Let hi{x) = 

When d{x) = we 


djx-v.) _ x-u 
d{v--u) v—u 

lave 


d(x-v) __ x-v 
d(u—v) u—v 


Note that 


|g7(w-x)/2 _ 1| \(,-l{x-u)/2 _ 1| < gl7|('u-ii)/2 _ g- 


|7|(t)— u)/2 


( 2 . 9 ) 


By applying the mean value theorem for the function on we got 

|g7N-x)/2 _ ^ |g-7(x-«)/2 _ l| < _ ^)|^| g|7|i/2. (2.10) 

Therefore choose C4 = I7I in this case. 

Suppose d{x) = e“'>'®/^ sin(/ 5 a;)// 5 . Note that in this case z = 1^1 Therefore 

-7(x-.)/ 2 sin (^(.7; - tQ) _ x-u 
sin (/?(v) -- •(/.)) 'I! — 'll, 

I sin(/:l(. 7 ; — ?/,)) x — u 


hi (. 7 ;) 




sin(/ 3 (u — u)) 


V — u 


fl !^lp- 7 (a:-w )/2 


( 2 . 11 ) 


H |e 

V — u 


11 - 
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Since is a strictly decreasing function in (0, z) we have 


sin(/ 5 (rr — u)) ^ x — u 


sin(/5(t) — u)) V — u' 

Therefore 

/..(x) < I - 1|. 

sin(p(n — uj) 

Similarly we have 

\smWu-v)) u-vj^'^ ‘I- 

Thus by ( 2 . 10 ) we get 

h(x) < - lU hi eW'/^ 

[sin(p(u — u)) sin(p('u — u)) J 


By Lemma 2.4 wo liavo 


COs(/?(7J — u)/ 2 ) 

< pj7l0-.0/2 [ 2sin'^(^(n-7x)/4) ) 


1 1 + {v — u) I7I 

l7|(/2 


Choose C4 = (^8 ^o^(Jl/Yj + lo'l) 

Similarly using the fact that sinh(a:) < xe‘, one can prove that there exists a 
constant C4 such that h(x) < Gtlv — u) in case d(x) = 2^;^^^ sinh((ri — r2)x(2). 

In case 7 = 0 we have 

h(T) - oSinh'((ri - r2)(u - u)/A) 

for all possible cases of d{x). Therefore there exists C4 such that h{x) < ^(d — uy. 
This proves the leimna. □ 


2.3 Subdivision Algorithm and its Refinement 
Equation 

We define the GBn-eurvt! a,s the Function 

n 

F{x)=^Y^Pil 3 nx;a,b), xe[a,b] 

i=:Q 
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where Pj € d > 1. To generate the curve F it is enough to generate each of 
its components. Thus we now present the subdivision algorithm for the evaluation of 
GBB-polynomial 

n 

fM = J2q‘lnnx-.a,b), .T£|n.i.l (2.12) 

1=0 

where G R, i = 0, 1, . . . , n and 0 < b — a < z. 

We denote J^nd by bo{l) and bi{l) respectively. 

Algorithm 2.1. 

Given the. control ponih qf, i = 0,1,..., n, loe obtain the control points qj, i = 
0, 1 , . . . , 2?! 05 follows: 

Let P’1 = qf, i = 0, l,...,n. 

For k = 1, . . . ,n, set 

Pi Pi~^^ i = k,k + l,...n. (2.13) 


Define 


Pit ^ 0 , 1 , . . . , Tl , 

i = n,n+ 1,. ..,2n. 


(2,14) 


In particular, p] — 6o(/) + b\{l) p® and 


pf bfil) pI_i +bi{l)pl 

~ (^'(0)^ Pi~2 + 2 bo{l) bi{l) + {biil))"^ p^i- 


By inducLiou wo got 


k 



j=0 






A: = 0, 1, . . . , n, 


k,k + \, . . . ,n. 


Therefore the set of control points qj, i = 0,1 . . . ,2n defined by (2.14) can be written 
explicitly as follows: 


<11 = 12 ( J ('^o(0)"'^ Ohior, k = o,i. 


, n 


s=0 

n 


'n - k 


(2.15) 


<iUk = Y2’i" („_ J (^'l(or^ 


Sz=:k 


A: = 0, 1, . . . ,7Z. 
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After the 7'tli iteration, the control points {(/fc'''’, /s = 0, 1, . . . ,2'^'^^ n} are obtained 
from {gjj, A: = 0, 1, . . . , 2’' n} as follows: For each j = 0,1, ... ,2^ replace ql by qjn+k> ^ — 
0, 1, . . . , n in the above algorithm. Then the resulting points in (2.14) are ^ ~ 

0, 1, . . . , 2n. Subsequently, we get 


= E (*) ('>o(4))‘-‘ ((>i(4))*, s = 0, 1, . . . . 

5 = 0 


n 


n — k 


I 


I 


(2.16) 


“ E 1" _ 'I ) » = n, 1 n. 

s=zk \ *- / 

In case I = b — a, the algorithm satisfies the following refinement equation: 




Yli=Q ^) ^)) X G [a, ^], 

. Z1"=0 9n+t ■^i‘(^; ^)) ^ ^ [^) 

which is immediate from the following theorem. 


(2.17) 


i=0 


Theorem 2.8. Let a, b E M. he such that 0<l = h — a<z and x' = Then for 
€ R, f = 0, 1, . . . , n we have 


^p’ n’‘(T.-,a,h) = 


Er=oPi a: G [a, a:'], 


i=() 


, ELo Pn~‘ ^-^"(a-'; a;', b), X G /;], 
where p^, k = 1,2, . . . , iv, i = k, k + 1, . . . ,n are defined by (2. IS). 


.. (2.18) 


Proof: We inovo the theorem by induction on n. When n = 1, 


d{x — b) 0 d(x — a) 


£p»B‘fca,6)=pS^^+p! 


Since 


d{b — a) ' 

u ,,a — b.d(x — a) ,, ,.d(x — x') 

d(. - « = rf(— ) + d(a - ^ 


for X G [a, x'] we get 


0 Di/ Po f ,a-b d{x-a) d{x - x') 


(2.19) 


( 2 . 20 ) 


+ Pl 


od(x - g) 


d{a — x') } ^ d{b — a) 


'^(V) ..0 , d(x-a) Qd(x-x') 


d(a-b) ^'“■^d(6-a) '''I d{!^) 


Pi 


TT + Po 


d(^) 
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d{x — x') 1 d{x — a) 


= Po-^J7^^+Pi 


i(<^) ' d(i^) 

1 

Y^p\ Bf(x;a,x') 


2=0 


In a similar way by expressing d{x — a) in terms of and 2(^Zb ) ’ S®*" 

1 1 

y~^Pi Bl{x\a,b) = B^{x-,x',b), x € [x',b]. 

2 = 0 2=0 

Tills proves (2.18) for the case n = 1. 

Let us assume that (2.18) holds for n = m — 1. We prove that (2.18) holds for the 
case n = m. By recurrence relation of i?-polynomials [38] we have for x G [a, x'] 


Y,p° Drix-.n.b) = ‘f^Pt Br'{x-.a,b) + B’r'(x-.a,b). 


'//I - 1 


m— 1 


_ g p; Br\x-a,x’) + Br'(x,a,x') 


d(-l) 

By induction assumption we get 

rri ./ 

53p“Br(i;a.6) = !i--- 
By (2.20) we get 

771 

Vp” nr(x-,a,h) 


d(l) 


2=0 


m— 1 


d(0 


2=0 


2 = 0 




(x; a, x') 


d{x — a) d.{-i^) . 1 ^— \ ^ i\ 

(x-,a,x) 


+ 


C?(^) d(b~a)^ 


x=0 


Rearranging the above equation we get 


III 

Y^plnnx-.a,!,) 


d{x - 


T7i- 1 


2=0 


‘^(t) i=o 


^pi Rf '{X]a,x') 


d(x - a) d{—^) . d{^) . „m_i, ,s 

d(0 ' (x;a,x) 


1/ Tfl 1 1/ \ TTl”! 

Br'(x-.a.x') + ^pSi Br'{x-.a,x-) 

2 1 t=0 

rn 


2-^ 2=0 


2=0 
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Similarly, by expressing d{x - a) in terms of and we get 

m m 

Br{x]a,b) = Br{x-,x',b), x € [x',b]. 

2=0 2=0 

Hence the theorem follows by induction. □ 

Remark 2.9. The above relinoment equation can also be proved using polar forms of 
a GBB function (see Remark 4.7(c) of [38]). 


After k iterations we have the following set of equations: for r = 0, 1, . . . , 2^^ - 1 

2=0 

where = a + r^. 


i=0 


Note that the function 


5fc,r(^) • — y ^ Qrn+i ^r+l)) ^ ^ [®r) ®r+l] 


( 2 . 22 ) 


2=0 


(2.23) 


is itself a GBB function defined on and has the control points 

where h k 

= 4 + r = 0,l,...,2^-l,.i = 0,l,...,n. 

77 ' 

For 0 < " ^jrl < n’ define flf{x) be the function in £„ which interpolates 

and at and respectively i.e., 

fT^kf \ ~ ^J,r)) k I d(n(x’ — K 

ci(n(e- - ej;.)) - e)) 

Let f^{x), X e [a, 6] be the function such that 


(2.24) 


^x) 






(2.2.'i) 


Then the curve {(x, /‘‘{x)) : x 6 K', ar+il) is the control curve for the GBB hinct 
Sk,r- Both the curves {(x, f^{x)) : x € [a]?, and {(x, Sfc,r(a::)) : x G K, 

lie in the D-convex hull of degree n of the corresponding control points 

{(4>> (Irn+i) ■ i = 0,...,n}. 
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Moreover SkM) = Qrn and Sk,Mr^i) = 9fr+i)n- ^ 

In Section 2.6 we show that the sequence of functions {f^} converges to the function 
f{x) on [a,i]. 


2.4 D-Convex Hull 


Let to < ti <•■■< tm, tm — to < bi € R, i — 0,1, ... ,m and 

D := {{ti,bi) : i = 0,1,. . . ,m}. 

In this section we show that the D-convex hull of the set B is the region enclosed by a 
simple closed curve. This curve consists of graplis of functions from £„ interpolating 
some elements of B. 


For any two points {U, bi), (tj, bj) of B we define an £„-function by 




djnjx - tj)) d{n{x - U)) 
d{n{ti - tj)) ‘ d{n{tj - U)) 


X 6 [to, t,fji 


and the corresiionding £„-curv(!: 


(2.26) 


• — {(•*') 9i,ji^'')) • ^ 


(2.27) 


Note that gij e C,^ and interpolates the points {ti,bi) and {tj,bj). Below we always 
assume i < j in the (Udiiiition of //jj. 


Definition 2.5 (Edge curves). Let 0 < r < s < m. A curve Gr,s is called an upper 
edge curve (UECj of B if 

fJr,a{U) ^ G [£)^a] 

and Ov.sifi) ^ '^ti ^ [tr? tj] , 

We call Gr,ii a lower edge curve (LEC) of B if 


grAk)<bi Vti€[lr,Q 

and 9r,s(f'i) ^ l^i '^ti ^ 


(2.29) 


We show in this section that the P-convex hull of B is the region enclosed by upper 
edge and lower edge curves. 



Chaptej' 2 


26 


Lemma 2.10. There exist contimious curves U and L defined on [to,tm] which are 
union of UECs and LECs respectively. Moreover union of U and L is a simple closed 
curve. 

Proof: First of all we observe that there exists U] e {1,2,..., m) such that Go.-ui is an 
UEC. 



^0 tp 

Figure 2.2: Existence of upper and lower edge curves 

Indeed, if the function satisfies <7o,i(4) > i = 2, 3, . . . , rn then (?o,i is an UEC. 
Otherwise, there is a least p € (2,3, . . . ,m) (see Figure 2.2) such that 6p > 

if Qo,i{^p) — K Hf'-iiiark 2.3 £fo,i = 9a, p and hence 

9a,p{U) > K i € {2, 3, ... ,71 - 1} and go.piU) = hi, i = 0, l,p. 

On the other hand if po,i (tp) < bp then po.i and po,p are distinct functions. Hence by 
Remark 2.3 

9o,p{ti) > 9i),i{U) > bi, ie {l,2,...,p- 1}. 

In any case we get po,p(^i) > bi, i e (0,1,..., p}. Now if go,p{U) > k for all i = 
p + 1, . . . , m, Go,p is an UEC. Else ■ there is a least pi G {p + 1 , . . . , m} such that 
.go,p{lpi) < bp,. Proceeding as before, we have 

go, Pi {li) ^ bi, i G {0, 1 , . . . , pi}. 

Therefore in finite steps, we obtain some ui < m such that Go.m is an UEC. 

In case u\ = rn, U — Go^m- On the other hand if ui < m, proceeding as before and 
starting from /ti instead of at 0 we find 1x2 G {'Ui + 1, . . . ,m} such that 


9ui,ufiU)>bi, X G {iii,...,W 2 } and 9ui,u2{U) > k, x G { 1 x 2 + 1, . . . ,?n}. 


Chapter 2 


27 


As Go,ui is an UEC, we have (7o,ui(^u2) > ^ 1 x 2 = ^ui,u 2 (^u 2 )' However in view of 
Remark 2.3, Quumiti) > go.uiiU) > bi, i e {0, 1, . . . , «i}. This implies that Gui,u 2 
is an UEC. 

Now if U 2 = m, U = Go,ui U Gui,u 2 - Otherwise we repeat the above procedure to 
find the upper edge curve This procedure can be repeated to obtain in finite 

steps the UECs Go.^i, which constitute If, 

C = Go,ui U Gui,U2 U . U Gur,m- 
An analogous argument establishes that 

L = Go,[i U Gj,,i2 U . . . U Gi^^m 
where Go,u, ...,Gi^^rn are LECs. 

It remains to show that 17 U L is a simple closed curve. 

Clearly U{xq) = L{xo) and U{xjn) = L{xm)- Let xq G [tp,tp+i], 0 < p < m — 1 
and Xo 7^ Then there exist an UEC Gui,ui+i and a LEC Gi.^i._^^ such that 

2:0 e and xq g Since is an UEC and is a LEC we 

have (7ui,7ii+i (^p) ^ .t^’XiiJXi+i (^ 7 >+i) — and (tp+i) ^ l^p+ii 

and hence by Remark 2.3 we get g[j^ij^^(xo) < puj.uj+i ( 2 ^ 0 )- This proves the lemma. □ 


For X G [t(),tm]i i'^t '»■(■>') and l{x) bo such that (.7;,u(.r)) G U and (.x, /(x)) G L 
respectively. Let us define 

C := {(x,?/) : l{x) < y < u(x), x G [to.tm]}- 

Theorem 2.11. C ix the snidMedt V -convex set, of order n contnining B. 


Proof: First we show that C is a P-convex set of order n containing B. ' 

For (xo,yo), {xi,yi) 6 C define p(x): 

, , d(n(x — xi)) d(n(x — xq)) r. . 1 

~ d(n(xo-Xi)) d{n{xr - xo)) ^ ^ 

Then we need to show that /(x) < p(x) < u(x) for all x G [xo,Xi]. We show that 
p(x) < u(x) and note tliat l(x) < p(x) can be proved analogously. 



Chapter 2 


28 



Figure 2.3: The case p(^) > u(^) 

I 

Suppose on the contrary there exists ^ 6 (a-’o.a;!) such that p(^) > u(^) and hence 
there exist Xq E (xq,^) and rr'i 6 (C,Xi) such that p(xo) = ^( 2 : 0 ) and p(x'i) = u(x'i) 
(see Figure 2.3). Let us assume ( 2 ;'o,ii(rc[))) G (^,u(C)) € Guj,uj+i and 

( 2 :i,«(a;i)) € Gu^.u^+i for some i,J and k. 

In case Xg, x[ G [tuj,tuj+i], then by Remark 2.3, p = guj,uj+i, which contradicts 
P(0 > w(0- 

Now suppose x\ > tuj+r Since Guj,uj.i.t is an UEC, by using Remark 2.3 we get 
guj,uj+j(x) > Pufc,ufc+i(-'c) for all X G if tu, i„j+i and for all x G if 

Therefore = P(^i)- Besides p(^) > 9uj,uj+i{0 

and hence p and guj,uj+i intersect in the interval ( 1 ^, Xj). Similarly we can show that 
Pui,ui+, (a-’o) = p(-^'o) < guj,uj+i (■A))- Therefore p and guj,uj+i have a point of intersection 
in [.'rg,^). Hence p and guj,uj+i have two i)oints of intersection in [.Xo, Xi]. This leads to 
a contradiction. 



Figure 2.4: A typical P-convex Inill. Filled squares represent the control points 


Now to see that C is the smallest I>-convex set of order n, let us take any P-convex 
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set Cl of order n containing B and sliow that every element of C is also an element of 
Cl- Indeed, let (;co,2/o) ^ C, then the curve 


g{x) := 


d{n{x ~ xo)) d{n{x - to)) 
Do -r 


2/0 ) ^ € [^ 0 ) 


d{n{tQ - xq)) d(n{xo - to)) 
intersect at least one of the edge curves U or L, say at {xi,yi) which belongs to Ci 
(see Figure 2.4). Since Ci is a P-convex set, (xofijo) 6 Ci as well. This completes the 
proof of the theorem. □ 


Define for r = 0, 1, . . . , 2*^ — 1, 

Dk{r) := max , flj {^)) - {y,&j))\\ (2-30) 

. 2/ e [^u,r,^v,r] 

0<u<v<n, ^ , 

where flj is defined in (2.24). Note that for a given n and for any k > log 2 (^) 
the inequality f is automatically satisfied for r = 0, 1, . . . , 2*^ - 1 and 

0 < i < j < n. The following theorem is an immediate consequence of Tlieorem 2.11. 

Theorem 2.12. When k > log2(n2/z) and r = 0, 1, . . . , 2*^ — 1 the diameter of the 
'D-convex hull of order n of the set 

9m+i) : i = 0,l,...,n} 

is equal to Dk{r). 

2.5 Convergence of Control Arcs and Control 
Polygons 

In this section we study the convergence of {Dk{r)} and show as a consequence the 
convergence of the control polygons and control arcs. 

First of all, we show that for any e > 0, there exists an N such that Dk{r) < t for 
k > N and for r = 0, 1, . . . , 2*^ — 1. Towards that end we denote 

Q{k) := max |gf| (2.31) 

0<i<2^n 
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- q^\ (2.32) 

0<i<2*7i-l 

and obtain some estimates for Q{k) as well as S{k). 

Note that by Theorem 2.12, the D-convex hull II of order n of {(^°, <j'°)}5bo is a 
bounded subset of . Since all the A:-th level control points 

{(et, <Jrn+i) : t = 0, 1, . . . , n. r = 0, 1, . . . , 2*= - 1} 


and their intermediate! poitits lie in //, there exists a constant C(j independent of i and 
j, and the level k such that 

l/'i'l < c. 


and therefore, Q{k) < C(j\/k. 
For 0 < .T < z define 


T]{x) := { 


“ rf(-i) ^ <i(x) - 

• r d{-x/2) . d { x / 2) 1 1 d{x/2) ^ d(-x/2) 

“ d(-x) ^ d(x) ^ anu ^ , 

jf d(-x/2) , d{x/2) , . dix/2) d{-x/2) 

“ d{-x) ^ d(x) ^ aim . 

Note that t){x) > 1 and < r]{x) for all x G [0, z). 

For a fixed I, 0 < I < z, lot us define <7(0) = 1 and for j > 1 


1 , 

^ d{x) > 
d{-x) > 


(2.33) 


(2.34) 


^(i) ~ (2-35) 

Observe that a(j + 1) = cr{j)r]{^) \fj > 0. 

I 

Lemma 2.13. For a fixed I, 0 < I < z, there exists a constant such that 

cr(j) < Ca Vi > 0. 


Proof: Suppose < 1- Then r}{l) = 1. Moreover, by (hi) of Proposition 2.f 

we have ri{^) = 1. Therefore a(j) = 1, Vi > 0. 

In case 4^ + > 1 and !^ > ^ then ,(0 = 2%^. and by (ii) and 

(iii) of Proposition 2.6 we liave r/(^) = Y) > 1. Tlii.s implies tiiat 


a{j) = 2^ 


j d{^) 

d{l) 


Vi > 0. 
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Therefore by (iv) of Proposition 2.6 there exists a constant Ci such that 

(^U) < Cl Vj >0. 

Similarly we get a(j) < Cx Vj > 0 for the case > 1 and 

■ This proves the lemma. □ 


For 0 < x < / let us define 




Theorem 2.14. For k > I 


s(k)<{o{k)r-‘{c,Y:0^ + ^ 


(2.36) 


(2.37) 


Hence there exists a constant such that for all k 


for general 7 
in case 7 = 0 . 


Proof: We prove the theorem by induction on k. By (2.14) we have 


^(1) = „ max 1 = max {|p*+i - p 1|, |p" * - p" *1} 

0<i<2n—l 0<2<7i-”l 

By (2.13) for a fixed 0<z<n — Iwe have 

|„i+l -i| - I _„i| 

lPi+1 P.l I P,+ P.+ l P.l 

^A-i/ 2 ) , ,7 <i(i/2) < i 4 i, , IpU.-p!I 

- “ 2 ' + 2 

^ ,nn . IPi+l-Pjl 
< a(0 C, + , ^ 2 . 

For a fixed i and j = i,i + 1, . . . ,n — 1 we have 

iPj+l Pj\-[ '■ 

Since, ( ^ P(0 P(0 > ^ rejuiating backwards we g(!t 
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Therefore 




pI 


<C,a(l) + {,il)y 


s(o) 


and 

o<K^^ll - pT'I < Oga{l) + (77(0)”"^^^. 

Since a(l) = r}(l) > 1 we get, 

5{i)<wi)r-'{c,«(i)+^}. 

This proves (2.37) for A: = 1. Now, let us assume that (2.37) holds for k = rn. Therefore, 


• 5(m) < 0^ + M}, (2.38) 

i=o 

Recall that i = 0, 1 , . . . , 2"‘‘''^n} is obtained from {q^]i = 0, 1 , . . . , 2'^n} by the 

subdivision algorithm by re[)lacing I by 0. Therefore, repeating arguments used for 
proving (2.37) for A: ~ 1, we have 


S(m + 1) < + 5(™)/2} 


< (-7(5!;))"-*^, 0(5!;) + Mrn + l))"-‘ S ^ + (''("• + l)r‘^ 


2 z— ^ 2^-^- 

j=o 


Since, 1 < rj(^) < a(m + 1) we have, 


S(m + 1) < (cr(m + I))"" ^ {Cg ^ 

i=o 


«(g) 5(0). 


This completes the induction and hence (2.37) is true for all A. By Lemma 2.7 and the 
fact that cr(A:) < Co- we have. 


rr/i ^ ^ f ^ 2/m <S'(0)., - 

S(k) < C0Cg— + for any 7 


and 




This| jc j» jidetes the proof of the theorem. □ 


for 7 = 0. 
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Our aim is to show that for ovary c > 0 , there exists an N such that for every 
r G {0, 1, . . . , 2*^ — 1}, Di^{r) < e for k > N. Towards that end, we introduce the notion 
of A:th level control polygon as follows: For any two points {^Ir^Qrn+i) (^j>. 9rn+j)> 
we define 

tk ck 

'■y (-x) := <i..H + <7t.+,-, X € Kt.eM (2.39) 

Sj,r S) 


!»z,r 


CK tK 

^'Lr Si' 


and /i*(a:) by 


Then, the piecewise linear curve {(a;, /t*^(x)) : x G [( 2 ,i»]} is called the A:th level control 
polygon of /. Now we are in a position to derive a rate of convergence of {Dk(r)}. 

Theorem 2.15. For k > log 2 ( j) and r = 0 , 1 , . . . , 2 *^ — 1 , there exists a constant C 
such that 


Dk{r) < 


C_k 
2k • 


Moreoxier if 'y = 0 then there exist a constant C such that 


'Mr) < 2 


Therefore limjt-^oo Dk{r) = 0. 


Proof: Let k > log 2 ( 7 ) and 0 < r < 2* - 1. Suppose 0 < i < j < n and 
0 < u < 'w < n, and x G ‘>''1^ V ^ [^u,r. we have 

- (y. = 1^ - yl^ + mfi^) - /uivWP 


< 


p 

22k 


+ - my)\ 


Tykt 


(2.40) 


Moreover 


(2.41) 


< l/'f (^) - - KtMl + KM - 

Since hl'j{x) and fdfk^{y) lie in the convex hull bf points {^rn+JLo middle term 
becomes 

(^) - Kviy)\ ^ ” o<T<^n^-l ^ 

By Theorem 2.14 we get, 

nCs for any 7 , 


|/^IJW-/Cf.(y)l< 


“nC, when 7 = 0 . 


(2.42) 
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Further 

er,k 


^ I rf(n(a; - 

d{7i{x - 


— F^ 

^3,r Si,: 


d(n(x - ^ 

Qrn+i "f" "JTCTck _ ck \\ ^rn+ j ~ Jk _ (k ^rn+i ~ Tk _ ck ^rn+j | 

^?.r Si 




z,r 


<a 


^j,r - 3^ I I din{x - <ef,)) X - 


d{n{^lr - ^Ir)) d{n{^l^ - ^l^)) - ^Ir 

By Lemma 2.7 we get, 

|/'f(x)-AyW|<C,C,^. (2.43) 

By substituting the estimates (2.42) and (2.43) in (2.41) and by (2.40) we get 


\{xJl;-{x)) - {yJui{v))\\ < 


where C 


C^, for any 7 
C^, for 7 = 0 , 

P + nCs + 2Cq C 4 1. This proves the theorem. □ 


Remark: Thus it is clear from the previous theorem that the sequence {/*(a;)} con- 
verges to f{x) uniformly on [a,b]. Besides, in view of the inequality 




K'-ix) 




< CqC,-^ 


(see (2.43)) the sequence also converges to / uniformly. 


Note that /(3:)|[„j;,aA:^^] = Sk,r{x)- Therefore f{x) is a X>-polynomial spline function 
as discussed by Gonsor and Noamtu in [38]. Hence l)y Theorem 7.2 of [38] the following 
important consequence of quadratic convergence of control curves follows. 


Theorem 2.16. For k > i) and x € [a, 6 ] 


ir(:e - /(.x)i < K- 


p 


where 


rr 


K = -- sup 

rh,...,Vn^ [«, 
1 <i,j <n 


22k 

d'^F 


d7jidy^^ 




(2.44) 


(2.45) 


and F is the polar form of f [38]. Hence the sequence of control curves {/*} converges 
to f uniformly and quadratically in [a, b]. 



Chapter 2 


35 


2.6 Applications and Examples 


In this section we present some examples, illustrating how the above scheme can be 
applied to modify the shape of the curves and surfaces by suitable choices of d{x). This 
results in groat advantage for design purposes. 


In case of classical Bezier curve: 

n 

Y^PiBi{x][a,b]), xe[a,b] (2.46) 

1=0 


the tangent vectors of the curve are fixed and are n(Pi - Po)/{b — a) and 
n(Fn — Pn-i)/{b — a) at the end points Pq and Pn respectively. However, for a given 
d{x) the subdivision algorithm generates a GBB curve having tangent vectors: 


cv = n(Pi 


1 


_ d'{a — b), ^ d'(b — a) _ 

+ Po-7r ~)y ^ — '^{Pn— — + Pl 


d{b — a) d{a — b) 


d{b — a) d{a — b) 


) 


at the end points Pq and Pi respectively. Consequently, for different choices of d{x), 
this algorithm generates GOO curves with different slopes and curvatures near the end 
points. Therefore, appropriate choice of d{x) can be used for designing curves and 
tensor product surfaces. In this sense, one can consider the parameters of the equation 
(2.1) leading to d{x) as design parameters. We illustrate this with some examples. 


In the following examples we take a fixed set of control points. 


Po = (1, 1), P, = (0.5, 2), P 2 = (2.5, 2), Pa = (2, 1) 


and I = 0.5. Four iterations of the subdivision algorithm are used and the control 
polygon at fourth level is used to represent the GBB curve. Different GBB curves are 
displayed in Figure 2.5(a)-2.5(f). Corresponding d{x) and the tangent vectors at the 
end points of the GBB curves are listed in Tal)Ie 2.1. 

The subdivision scheme for curves is extended to the tensor product scheme by the 
usual procedure. One advantage of this scheme is that we now have different design 
variables d{x) in the diderent directions. 
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Figaro 2.5: GBB oia'voH witli different end point tangent vectors 


Table 2.1 


Figure 2.5 

Tangent at Pq 

Tangent at P3 

d{x) 

(a) 

(-3.6) 

(-3,-6) 

X 

(b) 

(-171.88,343.79) 

(-171.92,-343.79) 

sin(rc) 

(c) 

(-3.61,5.02) 

(-1.41,-5.02) 

sinh(.r) 

(cl) 

(-1.70,3.75) 

(-39.82,-13.53) 

xe^ 

(e) 

(-236.13,76.36) 

(-722.18,-785.62) 

e® sin(rE) 

(f) 

(-1.75,3.49) 

(-4.75, -9.5) 

sinh(a;) 


To ilhisirate tlic tensor product schoine we take a set of control points suitable for 
a tensor product scheme. They form vertices of a polyhedron which is shown in Figure 
2.6(a). Figure 2.6(b) represents the classical tensor jiroduct Bezier surface where as 
the Figure 2. 7-2. 8 sliow dilfercnit tensor i)roduct CJBB surfaces obtained by using the 
above subdivision sclienie. We take only three iterations along both X and Y. The 
values assigned to d{x) along A" and Y directions are shown in Table 2.2. 
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Table 2.2 


Figure no. 

d{x) along X-axis 

d{x) along T-axis 

2.6(b) 

X 

X 

2.7(a) 

X 

sin(a:) 

2.7(b) 

X 

sinh(a:) 

2.8(a) 

sin(a:) 

X 

2.8(b) 

xe^ 

sin(a;) 




(b) 


Figure 2.6: Origina] polyliedroii and the classical tensor product Bezier surface 



Figure 2.7: Tensor product GBB surfaces for (a;,sin(a:)) and (i, sinh(a:)) 
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(a) (b) 

Figure 2.8: Tensor product GBB surfaces for (sin(x),x) and (xe^,sin(x)), viewpoint 
(- 78 , 28 ) 


Chapter 3 


Trigonometric Spline Curves 


III this chapter a noiistationary subdivision scluiiiK! for trigonometric spline curves is 
introduced. The scheme generalizes the Lane-Riesenfeld scheme for the polynomial 
spline curves. The chapter is organized into five sections. In Section 3.1 trigonometric 
splines are introduced briefly. The subdivision scheme is introduced in Section 3.2. The 
subdivision algorithm is based on a convolution formula which is derived in Section 3.3. 
The associated refinement equation is proved in Section 3.4. Convergence of the scheme 
is shown in Section 3.5. Finally in Section 3.6 the scheme is used to reconstruct circles; 


3.1 Trigonometric Spline 


Trigonometric B-splines {T”(a;; 1) : j = 0, 1, 2, ... , m} of order n (with the mesh size 1) 
associated with the knot sequence 

A. := {ti = il : i = 0, 1, m + n} , m> n, 0 < Z < ^ (3.1) 

are defined in [48] by the following recurrence relation: 


, 1 1 if a; G [0,Z), 

I 0 otherwise, 


(3.2) 
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and for n > 1, 


TSir.l) = + 5!’* - i;i) 


1 


s{tn - ti) 


s((n - 1)1) 
where s(x) — sin(a:). Moreover define 


{sWT„"-‘(x;i) + s(t„ - x}T„"-\x - );1)} 


(3.3) 


Tj'ix-J) :=T^{x - jl;l), for j = l,2, (3.4) 

The trigonometric B-spline T'^{x\l) is supported on [tj, tj+n] and 

{T;\x-,l):j = 0,l,...,m} 

are linearly independent, on the interval 'f'm+\]- It is known [48] that on [tn-b tm+\] 
any trigonometric spline /(;r) of order n with knot sequence A has a unique represen- 
tation; 

rn 

fix) = Y,n i). K 6 K, j = 1, 2, . . . , m. (3.5) 


Control curve and convex hull of trigonometric spline: 

For n > 1, let 

L„ := spaii{Hin(.(7), — cos((n — I).?;)}- (3-6) 

For the knot sequence A, let the moving knot averages are given by 


Then, the points 
metric spline /(.t), 


1 


i+n— 1 


t*: := — - — tj = (i + ^)l, 0 < i < m. 

* n - 1 ^ ^ 2 ^ “ “ 


(3.7) 


(t*, Pi), i — 0,1, . . . ,m are called the control points of the trigono- 
G- im-ti], defined in (3.5). 


It is easily computed that for any two successive control points (tj, Pi) and Pi+i), 
i = 0, 1, . . . , m — 1, the unique function gi{t) in 'hn interpolating pi and pi+j at t* and 
is given as follows, 


s((n- l)(il+i -i)) „ , 

■“ s((n - 1)((,V. - n)) s({n - l)(t(+, - *;)) 
The function g{t) on [t^, t*„], where g is defined piecewise by 

•^1 [t: , ~ 9i(l ’) ) 'I 0,1,...,?TI 1 

is called the control curve of the trigonometric spl ne function f{x). 


(3.8) 
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Definition 3.1. [Trigonometrically convex set) [48]): 

A subset B of is called trigonometrically convex of order n, n >2 if for any two 
points (^ 1 , Cl) and (^ 2 ) ('■z) of B with 0 < ^2 - the curve 




s{[n-l)[^2-0) s((n-l)(^- 6 )) 

1)(C2-6)) ^ 5((«- 1)(6 -6)) ^ 

Definition 3.2. [Trigonometric convex hull) [48]; 

The trigonometric convex hull of order n of a subset B of 
metrically convex set of order n containing B. 


■ ^ € 


( 6 ) 6 )} 


C B, 


is the smallest trigono- 


Suppose S[t) is a trigonometric spline function on an interval I with {(s* , as 

its associated control points. Then it is shown in [48] that the spline curve {[t, S[t)) : 
t & 1} lies in the trigonometric convex hull of order . n of {(s*, Ci)}”^]/. In particular 
the curve {[xj[x)) : a: G [i^n-i , i^m-ti]}, where f[x) = in the 

trigonometric convex hull of order n of : j = 0, 1, . . . ,m}. 


3.2 Subdivision Algorithm 


In this section we present the subdivision algorithm for the evaluation of the trigono- 
metric sidino function f[x) given by (3.5). We denote by w[j,l). 


Algorithm 3.1. Let n > 1 and n = 2A; — 1 or n = 2k for some k > 1. Given the 
control points qf, i = Q, I, . . . ,m ive obtain the control points qj , i = n — 1, . . . , 2m + 1 
as follows: Set 


P\i 



i = 0,1, . . . ,m 
(Qt-i + Qi)^ i = l,2,...,m 


if n = 2k — 1 
if n = 2k, 


Pzi+i = * = 0,1 ,..., m [n — 2k orn = 2k - 1). 


( 3 . 9 ) 

( 3 . 10 ) 


For j = 2,3, . . .k, define 

'w[2j - 3, 0 {w[2j - 2, 0 pi:^ + pill -b w[2j - 2, /) pi"'}, 
for n = 2k — 1, and i = 2j — 2, 2j — 1,. . 2m -f- 1. 

w[2j - 2; 0 {t(;(2j - 1; 1) + iCl + - U 0 }> 

for n = 2k, and i = 2j — 1, 2j , . . . , 2m -f 1. 
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Set qj = Pi, i = n - 1, . . . ,2m + l. 


The equations (3.9) and (3.10) are considered as the doubling up step and (3.11) 
is considered as the averaging step of the above algorithm. In the next theorem we 
show that for each i, q} is a weighted sum of <7°, j = 0,1,2, ... ,m. 

Theorem 3.1. For each i = n — 1, . . . , 2m + 1 the following refinement relation holds: 


= ( 3 - 12 ) 

'where / 0 onig for j = 0 , 1 ,..., n. These eoefjieienls are defined by the recurrence 
relation 


«•«.,(/) = 1 -«,.,(/), = 


j — 0, 2, 
j = 1. 


(3.13) 


and for n > 2, 




5((n-2)f)s((n-l)0 
s{{72 — 2)1) s{{n — 1)/) 


i )'*'s((n- 1)0 '■‘■"‘"0 


(3.14) 




Proof: We prove the theorem by induction on n. In case n = 1, by the above algorithm 
we have = qf and = r/°. Therefore 0.0,1 (/) = 1 = n,\,\{l) and ai,i(/) = 0 for 
i ^ 0,1. In case n = 2, by the above algorithm we have <72^ = ■^('/I’-i +9?) and 
92 i+i Therefore (3.12) holds true and ai,2(0 = 0 for i ^ 0, 1,2. 

Let us assume that the theorem holds for n = rn — 2 where m = 2k ov 2k — 1. Then 
for a fixed i, we have 

pf"' = 5^«i-2j>-2(09?- (3-15) 

jf6Z 

By (3.11) we have 

pf = u;(m - 2, /){'u;(m - 1, 1) p^Zf + pfr/ + w{m - 1, 1) pf"‘}. 

Substituting (3.1 p in the above equation and simplifying we get 

Pi = 5^ a;(m - 2, 1) |w(m -1,1) ai-2-2j,m-2(0 
j€Z 

+ 1,771—2(0 w{lTl 1, /) 2j, 771—2(0 ^ • 
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Hence (3.12) is satisfied for n = m with 

= w{m - 2, l)^w{m - 1, /)ai_ 2 _ 2 j,m (0 + ai_2j-i,m(0 + - 1> l)ai- 2 j,m{l)]ql 

This proves (3.14) for n = m. It is clear that = 0, i = 0, 1, . . . ,m. This completes 
the proof of the theorem. □ 


Repeated use of the algorithm outlined above for the generation of new control 
points, gives us at the r + 1th level the control points 

{,[+' : i = (2''+' - l)(n - 1). . . . , 2'+'(m + 1) - 1} 

from {f/£ ; i = (2’’ ~ I)(n — 1), . . . , 2’’(m + 1) - 1} by tlie rule: 

= (3.16) 

jez 

Clearly, this subdivision algorithm is uniform non-stationary and binary [35]. The set 
of coefficients 

■ ;= 0, ],..., n} 

is the mask of the subdivision scheme at the r-th level. 

So far we have presented an algorithm for the generation of new control points 
using a mask a^''\ It will be shown later in Section 3.5 that the following associated 
refinement equation holds, 

m 2771+1 j 

^ e (3.17) 

i~() i:r:n-l 

where {T-''{x\ 0 : i n — 1 . . .2rn + 1} are trigonometric B-splines (with the mesh 
size 1/2) associated with the knot sequence 

Ai := [ti^] = : i = n — l,n, . . .,2m + n + 1}. (3.18) 

At the r-th level of refinement, we get 

m 2’-(m+l)-l 

Y,q°TT{x-.l)= Y, 

i=0 i=(2’--l)(n-l) 


(3.19) 
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where {T^{x\ : i — (2^ — l)(n — 1) . . .2'’(m + 1) — 1} are trigonometric B-splines 

(with mesh size ^ ) associated with the knot sequence: 

■’= {ii,r = i = (2’’ — l)(n — 1), . . . , 2^(m + 1) — 1 + n}. (3.20) 

Thus we have obtained a subdivision algorithm, a corresponding set of masks and 
refinement equations. Hence we have a uniform non-stationary binary subdivision 
scheme(BSS). 

We note here that generalization of Oslo algorithm to trigonometric splines with 
uniform knots [48] satisfy the same refinement equations (3.17) and (3.19), and has the 
same doubling up step as in Algoritlun 3.1. However, the prc.sent algorithm is distinct 
from the Oslo algorithm in view of its averaging process. 


3.3 Convolution Formula for Trigonometric B-spline 


In this section we derive a convolution formula for trigonometric B-splines which will 
be used to prove the refiiKunent ecpjation in the next section. An analogous formula 
for classical B-si)line can be found in [20, 21, 42, 71]. The Fourier transform of the 
function / G L^(R) is defined by 

/ OO 

/(.x) d.x, teR. 

■OO 

For any m,n € N such that n — m is even we denote x{m)x{m + 2) . . .x{n) by 


The following lemma gives an explicit expression for (To^(-; Z))1[y). 
Lemma 3.2. For n> 1 we have, 

, n—1 s((v — 7 *)-) 

(To"(-; l)r{ij) = 2”(n - 1)! K{n) Hj ■ • ^ V ^ 

where 

A'(n,) : 


V- 3 


1 , 


K{n-1) 

s((n-l)0’ 


n = 1 
n > 1. 


(3.21) 


(3.22) 
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Proof. We prove the lemma by induction on n. For the sake of convenience we employ 
the notations Tq and Tq{- — 1) for Tq(-; 1) and Tq{- — 1; 1) respectively. It is easy to derive 
that 

{TiTiv) = - S{,A). 

y ^ 

This shows that (3.21) is true for n = 1. 


For n > 1, from recurrence relation (3.3) we have 


OT+'Hb) = ^ {{■^TSny) + (s((.+. - .) To“(. -my)} . (3.23) 


Since 


and ((/(■ — 0)^(y) “ any y ^ L\ we have 


(7'o"+')'ty) 

_ j_rPoT(^-i)-{r„")~(y+i) , e-'i!'-”)' (To”)^ + 1) e-‘(m‘(jonvmi 

i(ni) I 


s{ni) 1^ 2i 

I f 1 — 


s{nl) 


2i 




+ 2 ! 
Q-i{y-n)l _ I 


2i 


2i 


(ToTCiz + i) , 




(3.24) 


By t aking n. = 1 wc get 


s{l) iij=-i s{y-j) 
which implies that (3.21) also holds for n = 2. 


Ktiny) = 2^ ^ R,.. 


Let us assume that (3.21) holds for n = m. Then by (3.24) we have 

s((y-y-l)|) 


1 


s{ml) 


2”*(m - 1)!^ ^^5 s((?/ + m)-) JT(7n) J][ 


■m+l ?/ — y — 1 


e e-i-(v+i)i s{{y - m)^-) 
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K{m) 

s{ml) 


— 1)! g-'("‘+i)y5 TT”^ i)^) 


n. 


y-j 


(y 


= K{m + 1) 2”^+^ ml fC ^ — 


y-j 


Which proves (3.21) for n = m + 1. Hence the lemma follows by induction. □ 


Now we present an expression for the trigonometric B-splines as a convolution 
product, analogous to corresponding result for the polynomial E-splines. 


Proposition 3.3. For n> 2, we have 




where, 


( n-2 s((n-l)a:) 

s{(7i-2)l s({n-l)l)> 


An(x;0 = < 


n—2 

s((n—2)l) 


s(n^l)(2l~a:)) 


0 < r < I 
l<x<2l 


[^ 0 , elsewhere. 


Moreover 

r;‘(x; 0 = A„(x; 0 * TT-^x; l),iGZ,xE R- 


(3.25) 


(3.26) 


(3.27) 


Proof. From the previous lemma we have 

{T^Tiv) ^ 2'^ (n - l)(n - 2) s((y - n -f l)|)s((y -1- n - 1)0 

(Vr^Tiv) ~ “ 2)^) “1)0 [y - « -f 1 ) [y-Fn - i) 

Let us denote 


hi{x) = X[o,;i (x) and h 2 {x) = X[o,il e 


where X[o,/] is the characteristic lunctiori of / — dearly 




g((!/-»-H) j) 

2 / — n-h 1 


and 


(Ihnv) = 2e 


-t(y|n-l)| 


.s((2/-|-yt-l)f) 
y -In — 1 


(3.28) 
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Hence 


(h,,h,r(y) = ? r*” + 


i/-fn — 1 y + n — 1 


Thus from (3.28) we get 

OTHy) = 


n — 1 n — 


s((n - 1)Q s((n — 2)1) 

Since T^(.; 1) and hi * /i 2 are continuous funcfciows weh.jive 


(iH'-hm (Trm 


Tt(x-, 1 ) = 


(n-l)( 7 i- 2 ) _2 


Oi iJii <b;"0(i;0 (3-29) 


s((n - 1)^) (.s(n— 2) 1) 

Let us now calculate [hi * h'i) (x) : 

nOO 

{hi*h 2 ){x)= / hi{ 2 l)hix—u)(Si 

J~O0 


/O 


l);x; 




[«,!) n[=t-< 1*1 


If 0 < re < 1 then [0,1] f| [x - l,x] = [0,a;|. TheEreforeswlien 0 < x < I, we have 


(/i, */i 2 ) 0 ;:) = dy 




= e 


— i(n —1) X 


-i(rt -1^ a: 


■g.i 2 Cii-' 1 )!/ 

i ! 2 (?t] — I) Jo 

^i 2 C«- l)c _ 1 
i 2 (ra- 1 ) 


_ sQi - l)ca:) 
rt - 1 ” 

If 1 < a; < 2/ then [0, l]f\[x — l,x]= [x - I, I ]. 'lker*ef()*re , when I < x < 21, we have 
{hi * h.,) (x) = e-*'’"-' > f •e'^ ("-'>» rf' 2 / 

Ji—i 


,i2 (n — 1)®; 


i 2(-ti— u) 


$—1 


_ g -•(»!- 1 )* 


(n — l)ir _ gi 2 (n-l){x— 2 ) 


t 2(m— 1) 


s((n-l)(2i — iD) 


n-L 
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If x € [0, 2/]^^ then [0, 1] H ^ = </*■ Therefore, when x 6 [0, 2/]'^, we have 

{hi * h2){x) = 0. 

This proves (3.25). The proof of Tp{x,l) = A„(x,0 * follows by a simple 

translation of the variable. □ 


3.4 Refinement Equation 


In this section we prove the refinernerit equation (3.17) associated with the subdivision 
scheme. We require the following lemma in which A„(x; 1) is expressed as a linear 
combination of its binary refinements A„(x; |), Aa{x — 5), A„(x — 1] |). 

Lemnici 3.4. For n > 2 and x 6 [0,2/], we have 




1)1) 


/ 


/, 


+ 


i'((n — 2 ) 1 ) I s((n 


l)i) 2' 2 ’ 2' 


s({n - 1)() 


;b}. 


(3.30) 


Proof: We observe that An(x;/), A„(x;|) e lUi so also their translations A„(x - 
il\l) and An(x - i^\ |). We also oljsou've that A„(x; 1) and An(x - 1; 1) are linearly 
independent on [0, /] while A„(.'i; - 1; 1) and A„(.x - /; 0 are linearly independent on 
[/, 21]. Therefore, ' 

A„(.r; /) = n„A„(.7;; ^) + a, A„(x - ^), .x G [0, /] (3.31) 

An(.x; /) = bc,An{x — ■x',x) + hAn{x — /; -), X G [/, 2/] (3.32) 

iL ui ^ 

for some ao, oi, bo and bi in R. 


By taking x = | in (3.31) we get 


n-2 s((n - 1)1) _ 71-2 s((n-l)l) 

s((n-2)/) ,s((n-l)/) “°s((n-2)0 s((n-l) 5 ) 


Therefore ao = 


s{{n-‘2) 2 ) s((n 1 ) 2 ) 
5((n-2)0 s{{n-\)l) 
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By taking x = I in (3.31) and (3.32) we get 

n — 2 _ n — 2 s((n— 1)|) 

s((n-2)l) ^^s((n- 2)^) s((n - 1)0 

n-2 _ n-2 s((n-l) 5 ) 

s{{n-2)l) - %((n-2)i) 5 ((n-l)|)- 

Therefore ai = = h- 

Simillarly by taking a: = 3^ in (3.32) we get bi = = Qq- This 

proves the lemma. □ 


We are now in a position to prove, the main result of the section, the validity of the 
refinement equation stable! in tlie following theorem. 


Theorem 3.5. Let n > 1 and n = 2k or 2k — 1 for some k > 1. For x € 
we have 


2m+l 


2 m+l 


Y, j)= E 

2 = 0 2 -- 71-1 2 = 71-1 


where p'- is defined by (3.1 J). 


[^n— 1 1 ^m+r 


(3.33) 


Proof. We prove the theorem by induction with respect to k. Since Tl{x\l) = 
(X-; 0 + T,' (. 7 ;; 0 for x G we haw; 


m • 771/1 J \ 2 m-fl t 

YViTl(.x-,l) = E«.“ (rifeO + Ti+ai;^)) = E P.’t'Oij). (3.34) 

t =0 i =0 ‘"^ 1=0 


Which proves (3.33) for the case n = 1, corresponding to /c = 1 for n = 2/;: - 1. 


For A; - 1 for the expression n = 2k, we have n = 2. Then our interval of con- 
sideration l)eeoines [/.|,/.„ii i] = [1, (m + l)^]. Since is also a trigono- 
metric spline of order 2 on with knot seciuence Aj, there exists cti G R, x = 

1,2,..., 2m + 1 such that 


2 m+l 


E j). 


X G [Ai,A„,+i]. 



Defi ng a: = (22 + 1)| and a: = (i + 1)/ = (2i 4- 2)| in the above equation we get 
A S' + <71°) '^2i+i = Qi respectively. This proves (3.33) for the case 

poir 


De 

Th 

ffl( 


it 

t 

t 


us assume that (3.33) holds for n = r where r = 2A; — 1 or 2k. Then we have 

2rnH- 1 


7n 


(3.35) 

1=0 i=r-l 

prove that (3.33) holds for n = r + 2 with r > 1. By (3.27) we have 


in rn 

<ii 1) = e 9? 0 • ')) = ') * E “< ') 


2 = 0 


2=0 


y induction hypothesis (3.35) gives 


2m+l 


m £.in~r J T 

E?? 'tT^ix-,1) = A.„(.T.,i) . E 4 5)' 


2-0 


2 =r— 1 


le expression (3.30) for Ar+.2{r:; 1), and taking Sj := tpo ^2 


q°T;*\x-J) 

( 1 11 I ^ 2 m-H I 

i Si Ar + aG ' f ; 2 ) + '‘>’2 Ar + 2 ( a ' - 2 ! 2 ^ ^ 2 ^ 

V ^ 2 =r— 1 

> m-fl r / / / /I 

E pf { a . 5 ) + *2 rrHx --;-)*, 3 "'+^(x - i ; -) 1 . 

= r-l ^ 

ir interval of consideration is [fr+i,tm+i] = [(^+1)G (m+1)^] and 5) = 0 

- 1)7, (m + 1)7] for i < r and i>2rn + 2, we have 

rn 

E«? ''r'(.'o'0 


2=0 


2Tri-f 1 


2771 


2 rri — 1 


'y Si Pi n ) + ^2 Pi -) + Si Pi Tf ^ 2 ^( a :; 2 ) 


i=r4*l 

2m+l 


2 -^ ' * I iTi N ' 2 ' 

i=r i=r-l 


E {.’lPf + .'<2lt,+.-lpf-2}'7T+*(93 2) 


i=r4*l 

2771+1 


= E 4"' h. 


i=r+l 


oves (3.33) forn = r + 2. Hence the theorem follows by induction. □ 
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3.5 Convergence Analysis 


III this section we discuss the convergence aspect of the subdivision scheme. The 
convergence of tlie subdivision scheme is proved using the convergence characterization 
of the subdivision scheme via the convergence of control curves (Theorem 1.1). An 
alternative proof is also given using the characteristic polynomials of the subdivision 
scheme, following a procedure in [.35]. 

We follow Definition 1,1 as the definition of convergence of a subdivision scheme. 
First we show the convergence of the subdivision scheme using the convergence of 
control curves. A control curve has been defined in Section 3.1 by (3.8). Let denote 
the control curve at the ?'th level of the refinement. For its construction we refer to 
(3.19). 


(3.36) 


(3.37) 


Clearly, the functions f),* ^ Z are compactly supported, continous and 
||V^i(-; ^)L ^ 1 ^ Mo’^cover, for x 6 it is easily seen that 

cos(^) 

spo A-* 











The spline function f{x) has the expression 

r(m+l)-l ^ 

Yj xe[tn-utmu] 

i={2’--l)(u-l) 

after the rth refinement. The points 

((•„</[), i = (2'-l)(n-l) 2'(m + l)-l 

with t-., = h±Ld j"th:hi- , L , r . are called control points and the function 

% ,7 f j, — i 

/’ (x) = ^ ^ (/i '*A(^) ^)) ^ ^ [^(2'’ — l)(n-l),ri ^2''(m+l)-l,r] 


i=(2''-l)(n-l) 


with 


f ^ ^ ^ [AL.i,r) A,r)) 




df) 

0, else 


is the corresponding control curve. 



Chapter S 


52 


and 


cos(— 


+ 



-4 1 uniformly as r oo. Therefore, 


^ uniformly in E. (3.39) 

iez 

Since i G Z are compactly supported, by the Remark in Page 51 of [29] the 

functions ^) , i e Z] satisfy the stability condition 


where h := G ^“(2’), (see (4.2) of [35]). 



(3.40) 


The following theorem is a consequence of Theorem 4.4 of [48], which shows that 
the sequence of control curves {/'"} converges uniformly and quadratically to / in 

[^n— 1) ^m+l]‘ 


Theorem 3.6. Let f and be defined by (3.36) and (3.37) respectively, then 

^ 2?^ccs(^) 

xvhere L := D'^ + (n — 1)^ and K is a constant depends upon n and f. 


Now we hav(^ all the necessary conditions satisfied for the application of Theorem 1.1 
of [35]. Therefore our subdivision scheme converges as a consequence of the theorem. 

We now discuss an alternative method for the proof of convergence of the subdivi- 
sion scheme using thci characteristic polynomials of the scheme. Every non-stationary 
uniform OSS is associated with as(!t of mask and a set of characteristic polynomials [35]. 

The characteristic polynomial at the rth level of our subdivision scheme correspond- 
ing to the mask {ai_„(”), i = 0, 1, . . . , n} is given by 

= ^ = e“. t = x/=T, leR, (3.42) 

j=0 

Clearly, 

a[’'\z) ^ i z (3,43) 

and 

= + 2,:(i) I + ^)(1 + a- 

The characteristic polynomials r4r^(2;), n > 2, are given by the following lemma. 


(3.44) 
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Lemma 3.7. For n > 2 




2~n4-l 

T" 


1 

. ~ n + l 


rijLi(l+e n^i(l + z)(l 4- 2 ), n = 2m+l 


n;T:T,7(-p^) I1;L, (1 + n - ^(1 + n = 2m ■ 

(3.45) 

where c{x) = co.‘i{x). 


Proof: By (3.14) 


c((n - 2)^ c((n-l) ^) 2'-+i ^ 

-(1 + c-”("-»)27Tr 2)(l + z) ai,’22('2)- 


c((n-2)2;^prf;((n- 1)5;^) 

By repeated use of the above recurrence relation with respect to n, we get 



2~n+l 

0"=/ 

2-"+' 

n”=i‘ 


n;=,(i + n," .( 1 + 

n" , (1 + e-'®-'>F!Fr + e«-»r4T 


n = 2m + 1 
n = 2m. □ 


By Example 1.1 if the characteristic polynomials at r-th level of a non-stationary 
BSS is 

bj[\z) = 2 JJ -(1 4- 62 ^ z), r e Z 4 ., Cj E C (3.46) 

j=i ^ 

then the BSS converges and the limit function is in C'''“^(R). Further the limit function 
is an exponential spline. 


We use the above in order to show the convergence of our scheme. For this we 
define 

-'i(2i - 1);, j = l,...,m 

, for n = 2m, 

i(2j-l)l, j = m + l,...,2m 
-2ijl, ; = l,...,m 

< 2ijl, i = m + l,...,2m, , for n = 2m + 1. 

^0, i = 2m 

Then the characteristic polynomials {aJ‘j(^;)} in our scheme has the form 

(3.47) 


rij-i <3r^) . 
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Dt'note by 1%. Then for all k we have 


Tl-1 


n— 1 


riA =0^“^ <(-!)! 


s{jl) 

Since for j = 1, . . . , n - 1, 0 < / < i, A: € Z+ the sequence 

{nr=o/^''} is a monotonically increasing sequence which is bounded above and hence 
converges. In fact 


ri- I 


i 




r=:0 


Therefore our subdivision scheme converges and the limit function is in 






Figure 3.1: Quadratic trigonometric spline: (a) initial control points. Limiting curves 
after three iterations: (b) 1=0.25, (c) 1=0.5, (d) 1=0.75 

Towards this end we illustrate the subdivision algorithm in Figures 3.1 and 3.2. 
Note that we now have variable I wliicli is now used ns a design variable. Taking 
different value of I we get limiting curves with different curvatures for the same set of 
control points and control polygon. 
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Figure 3.2: CiiV)ic trigonometric spline: 
1=0.25, (h) 1=0.5 

3.6 Generating Circles 



Limiting curves after three iterations (a) 


In CAGD wo sometimes reciuirc to construct circular parts of geometric models. In 
this section we impleiuent our subdivision scheme to reconstruct circles. For this we 
take n=3 and the vertices of a regular rrr-gon 

Vj := ^cos(i^), sin(i^)^ , j = 0, 1, . . . , m (3.48) 

as the initial control points. The characteristic polynomial now becomes 

= --r-rv-7TT(i + 

cos(2FFr) cos(^) (3,49) 

^ cos( ^ ) 2-^ _ + z). 

cos(^) cos 2 ( 2 rTr) 

By Example 1.2 if the points V}, j = 0, 1, . . . , m are taken as the initial control points, 
then the limit curve of the BSS associated with the characteristic polynomial 

^ - (1 + e~*^ 2 :)(l + e‘^z)(l + z), 1 = 2t/m 

is the circle with origin (0, 0) and radius cos(27r/m) inscribed in that regular m-gon 
whose vertices arc given by Vj, j = 0 , 1 , . . . , m. 


r A ^ 

*:->«) cos(^) cos(/) 


Note that 
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- 2-1012 - 2-1012 


(a) (b) 

Figm‘0 3.3: (a) The .siiuarci whoso vertices lie equidistantly on the unit circle (dotted 
curve) is considered as initial control polygon, (b) Limit curve after third iteration 



Figure 3.4: (a) The regular triangle whose vertices lie equidstantly on the unit circle 
(dotted curve) is taken as initial control polygon (b) Limit curve after third iteration 

Thcrofore wl.cn ( = tlic limit enrvo of tl.e BSS a.tsociated with the polynomial 
is circle with origin (0,0) and radius 1. Since the distance of origin (0,0) 
to Vj is 1 the resulting circle passes through the vertices V^-, y = 0, 1 , . . . , m. 
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Interpolatory Scheme 


Th(! (haphn' is organ izc'd ns follows. In Section 4.1 the subdivision scheme is intro- 
duced and the convergence! of the scheme is also studied. Some properties of the basic 
limit function are presented in Section 4 . 2 . In Section 4.3 we show that the scheme 
reconstructs a certain class of trigonometric polynomials. The quadratic order of ap- 
proximation of the limit function is shown in Section 4 . 4 . The scheme is used to 
interpolate a complex sequence in Section 4 . 5 . 


4.1 The Subdivision Algorithm and its Convergence 


Let us define a space of trigonometric polynomials T by 

T ;= span{l,c(a;),s(a;),c( 2 x),s( 2 x)}. ( 4 . 1 ) 

where c{x) = cos(o;a:) and 5(2:) = sin(o!x), 0 < a < tt. 

Let g e T and g(x) = ao + aic{x) -f a2s(x) -f a3c(2x) + a4s{2x). Then the sum 
03 -f 04 is called the amplitude of g. 


Suppose wo have a data set 

D = {(xo, f{xo)), {xuf{xi)), (2:2, f[x2)), (2:3, /{xs))}. 
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Then theie exist several functions in T interpolating the data set D. But it is known 
[53] that the function 

3 

L(j.) = f{xj)Lj{x) (4-2) 


where 


Lj{x) = c{ 


X — X. 


) n 

/c=0, kj^j 


,■ 4’-^) 


(4.3) 


is the unique function in T , which interpolates D and has the minimum amplitude 
among other interpolants from T. We call the function L{x) a Lagrange like inter- 
polant of the above data. 


Define Xj -- j, j = 0, 1,2, 3 and x' = 3/2 then 

1 sin^(Q;/4) 


Lo(.7;') = L^ix') = 

Ly{x') = L.,{x') = 


2 sin(Q;/2) sin(Q!)’ 
1 sin^(3a/4) 


Clearly, 


Let Wo = 


2 sin(Q:/2) sin(a:) ’ 
sin^(3a/4) sin^(a/4) , 


8in^(a/4) 

2 sin(a/2) sin(a) 


Hin((v/2) sin(o') Hin(fv/2) Hin(«) 
Then I 4- u;o = asinS^ls&a) 


L{x') = -woifixo) + f{xo)) + {^ + ^)o){f{xi) + fix 2 )). 

To define our nonstatioiiary scheme, for k > 0, we denote 

sin^(2T[^) 1 


Wk 


2 sin(^) sin( 2 ^) 16 cos2(^) cos(^) 


(4.4) 

(4.5) 

(4.6) 


(4.7) 


Some estimates of Wk which are useful in our scheme are given in the following 
lemma. 

Lemma 4.1. For k >0 and 0 < a < 7r/2 
(a) g > Wk > i 

/ij _ ±1 < ^ for some constant C independent of k. 


Proof: The inequality rok > ^ is obvious from (4.7). Observe that 

i < I 

16 (' 0 .s'‘^(,_jj, I g ) CO,s(,2j, | v) g (^COH^(2t-(;r ) T C()H(jipi;T)^ 


■Hlk 


1 

< - 
- 8 
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Below we establish the asymptotic equivalence of {Sk} and {5}. 

Theorem 4.3. The nonstationauj scheme {5^} is asymptotically equivalent to the 
stationary scheme {5}. Moreover, the limit function belongs to C'(R). 

Proof: We have 

(iCIj 

By Leinnia 4.1(b) we get |a5^i.^2^ - ax+ 2 fi\ < ^ and hence 

l|S‘-S||<^- (-l-ll) 

Hence “ *^lloo ^ schemes {S'*} and {S} are asymptotically equiv- 

alent. Since {.9} is a C’ stationary .scheme, by Theorem 1.2 the scheme {5*} is a stable 
and C” converg<mt scheme. Moreover, from (4.11) it is clear that 

CX3 

(4.12) 

Jt=0 

As the scheme asvsociated with S converges to a C limit function, by Theorem 1.3 the 
scheme {5*} also converges to a limit function. This proves the theorem. □ 


We illustrate the subdivision scheme in Figure 1. Initially we take a finite set of 
points (Figure 1(a)). The limit curves after three iterations of our algorithms are shown 
in Figure l(b)-l(d). 


4.2 Basic Limit Function 


The basic limit function of the scheme 


(iSfc) is the limit function of the scheme for tlie 


data 




1 i = 0, 
0 i/O. 


(4.13) 
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since cos 2 (.x) + cos(x) > 1 Vx 6 [0,7r/4]. This proves (a). Also note that 
Wir - — = —f ^ ~ 

16 16 V cos2(5^)cos(^) J 

^ (^2 + cos{^)^ sin^(2^) 2^2 ^ 

16 cos"'*(^j,) cos(jjp^) 16 cos2(n!/4) cos(q;/2) 2^''+'*’ 

Tlie lernina follows by choosing C -■ a^/{cos\a/i)cos{a/2)). □ 


Now, wo present the basic algorithm which is a nonstationary subdivision scheme. 


Algorithm 4 . 1 . Given the control points {p- G R, i = - 2 , - 1 , . . . , u + 2 } the control 
points i = —2, —1,0,..., 2^'^^n + 1} at level k + 1 are given by the following 


recursive relation: 

-1 < ^ < 2*^n + 1 

= -'^^kPi-i + (^ + Wk)Pi + (^ + '>^k)Pi+i - WkPi+2, -1 < * < 2*n. 


Remark 4 . 2 . If we take Wk = for all k, then this scheme coincides with Dyn 
and Levin’s four point subdivision scheme [ 32 ]. The point t;he 

Lagrange like interpolant of the points 


f.i — 1 I. s , i jt\/**hl k \/*"l'2 k M 
{( 2^: ( 2k ( ok ’l^i+2)}‘ 


2^= 


Note that the set of points at (A: + l)th level of the algorithm contains all the points 
at the kth level and some new points. Therefore the initial set of control points is 
contained at all levels of the algorithm. Therefore, the limit curve interpolates the set 
of initial points. 


Let us denote our non-stationary scheme by {Sk}. The mask of {Sk} at the kth 
level is a^^^ = , 03*^} where 


Ak) 


a_3 — 03 


(k) 


(k) (k) 

-Wk, «_2 = “2 


0, a<‘» = 1. all = af 


1/2 + tD/j, 


( 4 . 9 ) 


Note that the subdivision operator S, associated with the four point scheme of Dyn 
and Levin has the mask a {n .3, • • - , a;*} when; 


a-3 = <13 


-1/16, a_2 = 02 = 0, ao = 1, O-i = Oi = 9/16. 


(4.10) 
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Figure 4.1: (a) Initial data, (b) Limit curve for a = 0.25. (c) Limit curve 
for a = 0.5. (d) Limit curve for a = 1.0. 

By Theorem 4.3 the basic limit function belongs to In this section we derive 

some basic properties of the basic limit function denoted by F. Let 

D„ := {|;, j e Z}. (4.14) 

It is easy to check that restriction of F to Dn satisfies F(^) = p” Vy. 

Theorem 4.4. F is symmetric about the Y-axis. 


Proof: We prove this by induction on n. First of all F{-^) - F{-^) forn - 0 since 
_ p(^-j) = 0, j e Z by interpolation property. Let us assume that F{^) = 
F{-i^), ; € Z, A: = 1, 2, . . . , n. Therefore, F{^-) = F{-:^), V;. Moreover, 


2i I 1 
pizLl—) 

\ 2n-f l ^ 




.1) ( 1 /2 H- U!n)FU-) + (1/2 + Wn)F {'^-^ ) - WnF{'' ) 


2 “ 

-j + 1 


1 , \ ^ ^ —j 2 ^ 


■‘(VnF{—'-;~i — ) + (1/2 + ( 2,1 ) + ^( 2" 0 2" ^ 


ptJM±l). 

* v on I I ’ 


On I 1 
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Hence •^’(^) F{ for all j and n € Z. From the continuity of F we have 
F{x) = /' (--.t) V.X € R which completes the proof of the theorem. □ 

I h(’ translation b (• — A:), A G M of F is tlie limit function of the scheme for the 
initial data 

n 1. i = k, 

!'?=! (4.15) 

I 0, i ^ A. 




(a) (b) 

Figure 4.2; Basic limit functions: (a) for a = 0.25 and (b) for a = 0.5 respectively 




Figure 4.3: Basic limit functions: (a) for a — 1.0 and (b) for a — 2.5 respectively 


We show that F is a compactly supported function with support in [-3, 3] (See 
Figure 4. 2-4. 3). 
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Theorem 4.5. The basic lunit function F vanishes outside [-3,3]. 


Proof: Let us take to = 0 and define tn recursively by tn+i =tn + Observe that 
tn ^ -On ^ ^n-i- We claim that the restriction of F to Dn vanishes outside 

[~^n)^n]- Since F is symmetric it is enough to prove that F vanishes outside [0,t„]. 
We prove this by induction on n. 


Observe that from R.cunark 4.2 
2k 

2n+l'' ‘ ' 2^^+! 


2 A: , 2A: + 3, 

+ FTTt) = ^(^ri + ■ orTIi ) 


-Wn{F{tn + ^) + + ~^)) 

+ ( 1/2 + Wn) {F{tn + ^ ) + F{tn + ^ ^ )) 


(4.16) 


and 


/ 2A;+1. . A; + 2. 


(4.17) 


From the above equations it is easy to check that F{ti) = —WoF(0) ^ 0 and for A: > 0 

I 9 A' 

F{ti H — — ) = 0 and F{ti + — ) = 0. 


This proves our claim for n — 1. 


Let us assume F{-r) = d whenever x € D,, and x > t„. Then by (4.16) and 
(4.17) we get /'‘'(t„ 4 .i) = —WnF{tn) ^ 0 and for A; > 0 we have F(tn+i -h 2 ^) — ® 
F(tn+i + llTf) = 0. Since[-tn,^n] is contained in [-3,3] forallnandlimn-,ooHn,in] = 
[—3,3] we have the required result. □ 

Remark: Since F is compactly supported and continuous function on i?, it follows 
that F(- - k) are compactly supported continuous functions as well. If y{t) is the limit 
function intcu’polating the initial seqinince {?/(?).)}. Then clearly, 

A;+3 

y(^) = ^2/(;)F(t-i)= X) l/(i)i"(i-i) ' (4-18) 

jez 3=^-“^ 

where k is the greatest integer such that k <t. 
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4.3 Reconstruction of Functions 

I 

In this section we show that certain functions can be reconstructed by our scheme. It 
is easy to check that if pf = 1 for all i at fcth level then pj+^ = 1 for all j at A: + 1 th 
level. This shows that the hinction f{x) = 1 is reproduced by our scheme. Another 
simple consequence of this fact is that 

'^F{t-k) = l, teR. (4.19) 

jtez 

Therefore the translations of the basic limit function F form a partition of unity. 

The functions cos (ax) and sin (ax) can also be reconstructed by our scheme which 
follows from the following lemma. 

Lemma 4.6. Let k > Q and n > 0 be fixed integers. Let Pj = cos{j^), —2<j< 


2*^71 + 2. Then, we have. 



P 2 /' - 

and 

, 4 +' =cos((2i + l)^), -l<i<2'=n. 

Similarly, if Pj =sin(j^) then, 



and 

pW =sin((2. + l)^), -l<i<2‘n. 


Proof: We first prove the lemma for the first case: pf — cos(f^). Note that 
4+1 = p\ = cos('i^) = cos(2f— ) 

and 

Pw = -■'"fc (pti +Pi+2) + (1/2 + {Pi +Pi+i) • 

Since 

1 sin^(32i^) 

2 ”*" 2sin(2r^) sin(^) 

pf + p«.i =2 ™s(P* + i)5St)“»(5^> 


Pi-t +I'i+2 = 2 cos(( 2 i+ l) 2 ffj)cos( 32 j^.,), 


and 
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we get 


P'2i + 1 


_ / -2sin^(^) 2sin^(3^) ,2i^a 

V2sin(^) sin(^) “ 2^=+' 2 siii( 2 ^)sin(^) J 2 2*= 

_ 2 sin(^]^) sin(^) ^ 2i + 1 a . 

2 sin(2^) shi(.ji,-) ' 2 2^ 

= co4{2i+\) --■-). 

Analogously, = sin(^) and the expressions 


iti +Pi+2 = 2siii((2z + 1)^) cos(3^) 


and 


a 


a 


Ih -t- Ph i ■-= 2.sin((2i,+ 1)^) cos(^) 




lead to the following exi)ressions for /47+\’ 




2 siii( n) Miii( ' 2*' ' ' 2 siii(-^iq , ) sin(;jl, ) 


2sin2(32^) _ , a A . .2H^ a 




2 2*= 


) 


iV 


sin((2-/: + )• 


I'liis i)rov{‘H the lemma. □ 


Corollary 4.7. If the initial data lie on a graph of a function / G T and the values of 
/ are given on a set of equidistant points then the limit function of scheme (4.8) exactly 
reproduces the original function /. In particular if we choose a set of equidistant points 

O-TT 

pO = (cos(ik— ), sin(fc— )), A: = 0, 1, . . . , n 

on a circle, and a = 27 r/n, then the limit curve is the original unit circle. The Figure 4.4 
and 4.5, illustrate the above observation. 
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Figure 4.4: In (b) three equidistant points are taken as initial control points which lie 
on a circle sliown in (a). The control points and the control polygons after second and 
third iterations are shown in Figure (c) and (d) respectively. 




Figure 4.5: In (1>) foiir ecpiidistant points are taken as initial control points which lie 
on a circle shown in (a). The control points and the control polygons after second and 
third iterations are shown in Figure (c) and (d) respectively. 


4.4 Order of Approximation 


Finally, we present a result on tlie order of approximation of the above interpolation 
scheme. To state tlie basic approximation theorem, we assume that for a fixed n G N 
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let h — - and I — [ 2/j, 1 + /j]. Suppose that ^ is a C^(/) function defined on I and 

-2 < i < n + 2. 

The basic approximation theorem states that the limit function obtained by the above 
nonstationary scheme with the data g{ih), —2<i<n + 2 approximates g{t) with an 
error of 0(/i^). Now, we state the theorem precisely. 

Theorem 4.8. Letp^ = g{ih), —2<i< n+2, h = 1/n. Let f be the limit function of 
the nonstationary interpolatory scheme ( 4 . 8 ). If g e C^il) then there exists a constant 
C such that 

11 / - »L.|o,| := mg l/W - 5(^)1 < 


Proof: Let us define the (X)mpactty sui)portod functions Di, i eIj by 

Bi(i) := F{j^ - i), ie Z, 

where F is the l)asic limit function of the nonstationary scheme (4.8). By (4.19) 


5;]Bi(x) = l, Vis/. (4.21) 

i& 

Let us denote the inl.i'rval [/•//,,.///,] l/y Then for x € Ik,k+i hy,(4.18) we have 

k Vi fc+3 

/(i) = = E s{ih)B,[x). (4.22) 




i=:k—2 


Let us define 


. sm{ai{k + Z)h-x)) 0 sm{a{x - {k - 2)h)) . 

^ 5R35ft) 

It is easy to check that Q € T and hence by Corollary 4.7 and (4.18) for x € Ik- 2 , k+‘i 
wa get 

j^Q(ih)Fi--i). (4.23) 


It has been shown in [ 48 ] that if (j € then there exists a constant Ci depending 

only upon g such tliat 

Ii^ 

•^IiOO,4_2,A:+3 ■" ^^COS(5a/2) 


(4.24) 
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Moreover by (4.23) we get 




f-Q 


A:+3 


X; oWD^ix)- Y, Qijh)Bj{x)\\ 


ooJk.k^i H z — / ^ Woo.Ik k+1 

j=k-2 


< 5 max |lB,(a:) 


JI<5-» 


^"■2<j<A:+3“ ^ ^ ••oO»^fc-2,A:-h3 11'^ ^ lloO,/jt-.2,Jb+3‘ 

Since the basic limit fimctioius j € Z are bounded functions. We liave 

11-^ " - 11^ “ ^lloo./.,,+, + \\Q - i?IL, 

This complet.es tlie proof of the tlieorem. □ 


cos(5a;/i/2) ' 


(4.25) 


4.5 Complex Interpolation 


As an illustrat.ion of the our algorithm, we study the following problem of complex 
interpolation. b('(. c — re’*' Ix' a coiniilex number different from zero. Let us define a 
complex seciuence {^(/t)} by z{n.) — z"" Vn. Let x{n) = Re(z(n)) and y(n) = Im(z(n)) 
Vn. Let x(t) and y(t) be the C‘ limit functions of the scheme {5^} interpolating the 
initial sequences {.T:(/i)} and {;(/(‘n)} respectively. Then the function z(t) = x(t) -h 
iy(t), /. € R is called the complex interpolation of the sequence {^(n)}. 


Theorem 4,9. Let z{t) he the annplex interpolation of the sequence {> 2 (n)}. Then z{t) 
has the Fotirier series expansion 


^’(0 = E 




and 


where c„ UZ-a 

/,>,(,) .f. Ul±:^t± + 1 + (1/2 + w,.i)z - Wk-iZ^)/2. 


(4.26) 


Proof: Let us consider the function f{t) — Then f{t + a) j.t+a gib(t+a) • Since 

^(t+Q;) is the interpolation of the sequence {2(710+0;)}, 2:(na+a) = 2 ;^"'*' ^ =r e 2 

we have z{t + o) = r"c*'''’ 2 (t). Thus 

r“e**’“2(t) 

f{t + a) - ^c,^iba^tQibt ~ 
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I 

This implies that f is a periodic. Since z'(t) is continuous, f'{t) is continuous and 
hence 


/(*) = E “ 


(4.27) 


/:=“0 


where cjt = ^ /t“/(y)e dy. 


Define Ujt = ^ . It is easy to check that cq = /(O) = z(0) = 1 and 

^ ff \ I r /(O) + /(f) 1- 

Co = - f{x)dx = hm ^ — -t ^ = hm a^. 

(y Jq ^~>oo 2^ k-^oo 


Now 


/(^^) + /( 25 t) H ^ /(^) + /(3|f ) + •■•• + /(^^a) 

= ^ _ 

r 2*= - 1 


OA: 1 /(.;l ) I- i ('‘hT) -I 


Since /(p.ft) 


= -~4n- when p is odd we got 

=-■ --«^a-i2((p -- 3)^) + (1/2 + w;a-i)2((p - 1)^) 

+ (1/2 + Wfc-i)^:((p + 1)^) - Wk-xzi{p + 3)^) 


ivnd 


/(!>§) = + i(l/2 + w.-.)/((p- 1)^) 

+ 7(1/2 + Wfc-i)/((p+ 1)S) - 7Wi/((p + 3)^) 


where 7 = r . Thereforca 

f{w)+ + 

'"^Ok 


= J_ (zrt ± /(--+) + l+±^/(n) + (1/2 + Wk,,)/(5^)7 - 

2*= \ 7-‘ 2^’“' 7 z z 


Wfc-i 

/y3 


/( 0 ) + 


7 

1/2 + Wfc-l r , « 


7 


/(^) + (1/2 -H ^^^fc-i)/(2^)7 - 


+ 

4 - 1 /^ + ^ Qi ) + ( 1/2 + Wk - i ) f { oi)y 

ry3 2^'“^ 7 2^ 


+ 


7 



Using periodicity of / we get 


(Ik 


[I±l .p + 1 ^ (1/2 + ^^)fc_l)7-^0Jt_lT^}- 

2 2 7 '^ 7 


Which implies 


Ofc 


= afc_i/?fc(2fc,o) — • • • — 


j=l 


Hence r,, = Ihuk-yooak = evaluate c, k > 1 coasider the 

function co oo 

s(,) = c-¥“/(() = E = E 


j=-oo 


J=^oo 


Note that ,,(t) = "■'-o n. = 6 + f *• Therdbre by tollowiag the 

san.o , , too, ..lure udm.! r„ tvo pd e, = .'l, = R” . Kb(4.)' Sitte t',„ = z,, the 


theorem follows. □ 


Chapter 5 


Design of Surfaces 


A noa-s(, at, ionary .subdivision .schoinc For gonorating surfacc.s from arbitrary topologies 
is introduced. The limit surface is a bi-quadratic tensor product trigonometric spline 
surface except at a .small number of points termed as extraordinary points. The chap- 
ter is organiz(Hl into nine .sections. In Section 5.1 the topology of a data set is defined. 
The .subdivision scheme for a tensoi' product, topology is introduced in Section 5.2. 
The subdivision scheme for an arbitrary topology is introduced in Section 5.3. The pa- 
rameterization and prolongation procedures associated with the scheme are discussed 
in Section 5.4 and Section 5.5 respectively. The associated subdivision matrix is ana- 
lyzed in Section 5.6. Convergence of the scheme and continuity of the tangent plane of 
the limit surface arc shown in Section 5.7 and Section 5.8 respectively. Finally, some 
illustrative examples are presented in Section 5.9. 

5.1 Topology of a Data Set 


For designing a surface initially we are given a finite set of points. Usually the points are 
considered as vertices of a polyhedron. The lines joining the control points constitute 
edges and faces of the polyhedron. 

Order of the face: The number of edges consistuting the face. 



0rd6r of tho vertex 1 ho iininbcr of the edges eniarifiting from the vertex. 

A polyhcciroii is said to liave a teiivscr product topology if all its faces and vertices 
have order 4. If the orcku' of the faces and vertices are arbitrary, the polyhedron is 
said to liave an arl)iirai‘y topology. A data set have Iho to[)ology of the correspoudiug 

polyhedron. 


5.2 Scheme for Tensor Product Topology 


Wo recall l.hat in Chapter 3 a trigonometric spline curve is generated by a binary 
subdivision hcIk'Iih' (liSS). Usually a liniU; set ol’ cuiitrul j)uints (yj” G R, i € Z) is 
given. Tlien the BSS for a quadratic trigonometric spline recursively generates new 
control points which at A:-th level arc obtained by the rule: 


+ ui\Pi 

P2i-H =wipti +^oPt^ 


where xvq 




8in(4/i/2^) * 


(5.1) 


The above BSS is generalized to a tensor product BSS as follows; Let 

= {p"' e i,j e Z} 

denote the initial set of control ])oints which has a tensor j^roduct topology. These 
are aslo called control points at the 0-th level. Then the tensor product BSS defines 
recursively the A:-tii level control points 

{p^■€E^^ieZ} 


following the rule: 


P2i,2i 

P2i+l,2j 

P2i,2j+1 


= Wi {wipi^lj^i + U^oPij-l} + '^0 {‘^lPi-l,j + ^0Pi,j } 

= Wx + Wip'^jlx} + Wo [woptlj + ^ 

= lUo {wiPillj-l + WoPtjh} + + ^oPij'} 

= Wq {Wo PiCiJ-1 "b Wx Pgjil } + '^1 {^0 Pi-1 J *b '*^1 Pi,j }' 


(5.2) 



wlK'l'd 


M = 



8in(|) 

0 ^ 

8in(2/i) 

flin(2/ij 

Bin(|) 

8in(i^) 

0 

8in(2/i) 

sm(2k) 

V 0 

mn(^) 

siu(2/i) 

j 

8in{2/i) / 


Tlie sui)erscrii)t t in (5.3) denotes the matrix transpose. 
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Figure 5.2: Old and new control points 


(5.4) 


After simplification we have 


Boo = 


P + + jE/2 + F' 

4 ~ 


(5.5) 


where 


and 


P = 

El = 

E.2 = 


1 Aqo + A()i 

cos2(/i/2)cos(/i) 2 

1 Aqo d- A 10 

cos^(/i/2) cos(/i) 2 


(5.6) 

(5.7) 


(5.8) 


1 Aqo + A()i + ^10 + 

^ ~ cos2(/i) cos2(/i/2) 4 


(5.9) 


The application of tensor product schemes always requires that the data points lie 
on a tensor product topology, consequently, the smoothness properties of limit surfaces 
are consequences of similar results for the univaiiate case (see Chapter 3) 






Figure 5.1: Original polyhedron and limit surfaces after first, 


second and third iteration 


The above Figure illustrates the scheme (5.2). 


Si..cc a n.inimu,n of three co.rtrol points are required to gel a quadrat.c tr.^ ^ 
...otric splirre curve, a uriuhnum of ..iue control points with a tensor “P 

is req.iinnl to get a l,i-<|nadralie tensor |.ro.luct tngoiio.netric sp me pa o i . ^ ■ 

such a set of nine control poi.rts (see 1-igr.re 5.2) an.l denote then, b, A.,, J - 0, , . 

Using the Rule (5.2) four new control poi.rts eorrespondhrg to ea^ f«e 

sixteen new control points are obtained after the first .terat.on^ Let us d 

1 • • threG row B-nd first thrcG column by tj? 

new set of control points lying .n the first three row 

, ^ Tliis set of new control points has a matrix relation 

i,j = 0, 1, 2 (see Figure 5.2). ihis set oi nts 

/Boo Boi Bo2\ 

= M 


Dio Bn Bi2 

\jB 20 B 21 B 22 J 


Aio An ^12 

^^20 -^21 -^22 J 


M‘ 


(5.3) 




Quite often the coiitiol points lie on an arbitrary topology. The purpose of present 
investigation is to g(>n<u'ali/-e the above tensor i)ro(luet BSS in order to generate a 
surface from the control points lying on an arbitrary topology. 


5.3 Subdivision Scheme for Arbitrary Topology 


Let be the initial control polyhedron with an arbitrary topology consisting of planar 
faces which are convex. The subdivision scheme defines the new control polyhedron 
by the following procedure. 

Let be a face of F'* with vertices Aj, A2 , . . . , (Figure 5.3). A set of new 
vertices associated witli the fa,c(; 7'’” and the old vertices AjS is obtained by tlie sub- 
division rule: 

R1 : (li = ^ + + (5.10) 


where 


P 


Aj p, + Aj+i 


cos^(h/2) 


and 


F 


2cos^{h/2) cos{h) 
Aj T . . . + An 


cos^(/t) cos2(/i/2) n 

where all the indices are taken modulo n. After simplification we obtain 


a-i = a 


(n, 1) Ax + /3{n, 1) A.-i + /?(n, 1) Ax+i + 7(?b 1) (^i+2 + • • • + Ai+„-2) (5.11) 


where 


7(n, 1) 


1 


4n cos‘'^(|) cos^(/i) 


(5.12) 


a 


(n, 1) = 7(n.l) + + 4cos^(|)' 


Once all the new vertices associated with all tile faces of are obtained the 
following connectivity rule is applied to join them. 



(A) Tor a face h of join the now vertices associated with the old vertices in the 
same order as the corresponding old vertices. 

(B) Consider any old vertex of and the new vertices associated to it. Join each 
pair of new vertices across the adjacent faces of 

(C) Consider two adjacent faces Fi and F 2 in Let Pi and P 2 be two new vertices 
iissociated wit,h F\ and Fi respectively and associated with the same vertex of 

Join I\ and P^. 

As a consequence new faces enclosed by new sets of edges are obtained. Thus a new 
polyhedron P' consisting of these new faces is generated (see Figure 5.4). We denote 
the whole process by an oi)crator Fi, i.e. 



Analogously, another polyhedron (Figure 5.5) is generated from P' using the 
subdivision rule R2, where the subdivision rule Rk, > 1 is same as the subdivision 
rule R 1 but for the replacement of a{n, 1), ^(n, 1 ) and 7 (n, 1) by a{n, k), /3(n, k) and 

7 ( 11 , k) rcsiicctively, where for A: > i 

1 1 

(-08^(21^) 

1 

/3(n, k) = t(», e + 


1 



Here 0 < /(. < |. In general, Uie fc-(,h (control i)C)lyhcdrou is generated from ' 
using th«' subdivision rule Rk. This process is exi)ressed by 



Figure 5.4: 7^'^, the control polyhedron after first iteration 



Figure 5.5: 7 ^, the control polyhedron after second iteration 

Our ai.u is to study the convergence of the sequence to a smooth surface. 

First of ali note tluat the connectivity rule ensures that a face of 7 of or er m gives 
rise to a new face of order m; an interior vertex of ^ of order m generates a new fa 
of order » sorrounding it and each interior edge of ^ gives rise to a 
Moreover, ali the interior vertices in will have order 4. Therefore, " 

of quadrilaterals inere«s, nnmher of non-quadrilaterals 
iteration 5jt, > 2. However, after each iteration the jzes o ese 
Besides, the quadrilaterals sorround the interior non-quadnlaterals and 
these quadrilaterals inereases as we progress with the iterations. 

r 5 2 the layers of quadrilaterals sorrounding the 

Following the discussion of Section 0.2, tne layeis 4 

° 11 ii< 3 P nf subdivision rule and connectivity rule 

non-quadrilaterals obtained by repeated use o . 



(■()nv(’ig(' to a t,(ni.sor j)rodu(:l, l)i(}iia(lral.ic .s])lin(! siirlacc, while Uie iion-ciuadrilatcrals 
shiink to points t.einuHl as extraordinary points [26], I'lici ahovi! surface is C* oxci'pt 


possilily at, a small nunilior ol extraordinary points. In the following sections, our main 
concern shall be the study of the continuity and smoothness of the above limit surfaces 
at these extraordinary points. 



To study the nature of the limit surface near an extraordinary point it is enough 
to take .F” as a typical union of an n-sided face Fq and two layers of quadrilaterals Uq 
surrounding (Figure 5.6). 


5.4 Parameterization 


In this section we outline a procedure for the pararnetrization of the tensor product 
spline surface surrounding an extraordinary iioint. 

Let us consider a polyhedron as described above, hirst of all, we index the 
vertices of F° as 

Ziy j, i ~ 0, 1, . . . , 'U — 1 ) j “ 1 1 2> • ■ ■ 1 

(see Figure 5.7). The indexes 0, i and j corresponds to the level, the block and jth 
entry of the block respectively. These vertices also include the vertices of Uo having 
a rectangular topology. Therefore, a tensor product surface is obtained as a limit of 
and is denoted by X. The parametric representation of X follows the normal 
tensor product procedure, the details of which is given below. 


*'(V0 



Figure 5.8: Pat, dies A'o (dark shaded region), Xq^l and Xo,r (light shaded regions) 
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Figure 5.9: Parainetrization of the patches Xi, Xi^i and Xi.R 
Let J = 3n and J — {1, 2, ■ ■ • , J)- 

Parameter apace: The parameter space fi. = »(0) x J. of X consists of J copies of 
the square m(0) provi.lcd with the uciglihorl.ood relations as given in (5.14) eno e 
by := m(m) x J. J copies of the square »(m) = P, ^1 x [0. ^1, provided with the 

same iieighbo.irliood relations as {5.14) except that h is replaced by 

Representation: A representation of Q.m is a map 

4 , : (w.i) -> ,^(>") €K\j€j,we w{m) 

from n. to K^ where the functions ^ are smooth iniective maps in « the 
property that pieces of boundary curves identified by the respective neighborhood 



relat ions have idenUcal images. 


Now, coming t.o l lu' actual construction of we consider a bivariate affine map 
/ of the form 

fix, y) = {ax + b){cy + d) (5.15) 

which niajw a (luadrilaUu'al in into a ([uadrilateral in R^. It is easy to see that 
there exist affine maps /j, fi^n, * = 0, 1, . . . , n — 1 of the form (5.15) satisfying the 
neighbourhood relations (see Figure 5.10) 

fi{h,t) = fi,RiO,t) 

/i, /.(■*>, h) = /i(.s, 0) 

/i.z-(0,/i) = /i.ir(0,0) (5.16) 

fi-i,n{s, h) = fi,L{s,0) 

h) = /(),/.(■■*, 0 ) 


7 



Figure 5.10: Representation /o, fo,L and fo.R 

Clearly, the patches Xi, Xi,L, Xi,R can now be parametrized by the above maps on 
/i(«K0)). kiX’>‘’W) ntspoc.tively. The functions <l>\ ie J are obtained 

by re-indexing the maps ki, fi and k^^ * ~ 0, l,...,n- 1 sequentially. Ihus, the 
map (see Figure 5.12) 

(p : {w,j) <i>^{w), j €: J , w €. ■a'(O) 

is a representation of the parameter space Ho- Therefore, X is a smooth function on a 
mesh Fo (see Figure 5,12) defined by 

r„ = u"^pA’0n(0)). 


(5.17) 



Figure 5.11: The surface A'' 



Figure 5.12: Parameter domain Fq 

5.5 Prolongation 


We have seen in the previous section that tends as A: — > oo to a surface X with 

a hole inside it. In this section we see by the method of prolongation that the surface 
limA;_+oo is a smooth surface without any hole and is obtained by prolonging the 
surface A'. 

Before considering the nature of the limit surface obtained by limjt_^oo let us 
consider the p(dyhedrou X' = Si[T% It couHists of an 7i-sided face Fi ami four layers 
of quadrilaterals (Figure 5.13). The polyhedron - Fi has a rectangular topology 
and contains Si[U°] which is already known to have a rectangular topology. Therefore 
the corresponding limit surface (Figure 5.14) 

A' = lira SP [A' - Fi] 

k-^oo 




Figure 5.13: Control polylicdroii 



Figure 5.14: Limit surface X' 

The surface liiiifc-^oo 5'^*^ is same as X, obtained by applying the subdivision 

scheme {Sk}k >2 to the subdivided polyhedron -Fiff/®]. 


The scheme littifc-^oo U^] is now explained below in detail. 

For the sake of simplicity we consider 17® as a polyhedron having 9 control points in 
a tensor product topology. A single patch X is obtained as lim^-^oo'S'i The patch 
X has parameters defined on w(0) (see Figure 5.15). 

The polyhedron SiU^ is again a polyhedron having 16 control points in tensor 
product topology. From SiU^, we take 9 control points each (shown inside the dotted 
boxes in Figure 5.16) in tensor product toi)ology and obtain the corresponding patches 
by the scheme limfc-^c» • These patches have jnirameteis defined on m(l). 




Figure 5.15: as A; -4 oo 



Figure 5.16: The scheme as k oo 


Because of the nature of tensor product schemes and subdivision schemes, 
limjfc_,oo is a surface which is union of 4 patches with a neighbourhood relation. 

This surface is same as the original patch X but their parameters are now defined on 
a domain which is union of four w{l) rectangles with the same neighbourhood relation 
as their corresponding patches (see Figure 5.17). 

We now return to our original C/®. Following the above discussion the parameter 
space for the surface linik-^oo SiU^ which is denoted by is now obtained by 

dividing each unit of flo hy four equal parts (see Figure 5.18). The neighbourhood 
relations are obvious for this parameter space. 

Let <f)i be a representation of fli. Define the mesh 1 1 (figure 5.19) by 

(5.18) 


Then using the parameterization procedure ol Section 5.4, we see that the surface X 





is also a smooth surface on Fi. 



Figure 5.17; 5,77° as k oo 
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Figure 5.18: Subdivision of the parameter space 



1 

Figure 5.19** Puru-incter domuin Fi 

The mesh F. (Figure 6.20) is regular in the sense that each iuterior v«tex of F. is 
shared by exactly four patches. We assume that the regular mesh F, can be prolped 
by a layer pr(F,) such that S(F„) := F. Upr(r,) is still regular (see Figure 5.20). 

Using the parameterization procedure of Section 6.4 we see that the ^ ^ 

° , , r V' _ Y which IS a smooth exteu- 

smooth function on 5(ro) surface; (-^ ) • ^ pio-p.p 5 14') 

1 Af Y to orfFnl (light shaded surface m Figure 5.14). 

Sion of X called a prolongation of A to pr[i o) V g 



Figure 5.20: Prolngation of Fq 


One can continue the prolongation procedure iteratively. In general pr'^{X) - 
pr ’ (pr*^- ’ ( A' ) ) , A; > 2. Let P denote 

P := U (5.19) 

Below we investigate the nature of this surface P. 

Dcflnitiou 5.1. [63] A subdivision procedure {S\-} is said to be convergent, if thorn is 
a p such that for any seciuence, € pr’'‘'(A'’) we have linVa-^ x,n = p. Therefore, 

if a subdivision scheme converges, then P = PU {/;} is a surface without gap. 

The point p described above is usually called an extraordinary point and the defini- 
tion of convergence stated above ensures that the surface P is necessarily continuous at 
p. The convergence of a subdivision scheme depends upon the nature of the subdivision 
matrices which we study in the next section. 


5.6 Subdivision Matrices 

I 

5.6.1 Matrix Formulation 

We transform the subdivision process into a sociucnce of matrix transformations, whose 
underlying matrices are called subdivision matrices. First of all note that each of the 
tensor product surfaces pr^(X), m > 1 is defined by 9n =: K control points. We 
denote the set of control points of by a 9n x 3 block matrix Bm consisting of 

n blocks. Here m represents mth level. 
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Figure 5.21: Control points at rn th and (rri + 1) th level of subdivision 
Using the indexing j)rocedure used for Brn is partitioned into n blocks 


^rn,l3 • • • j Bm, 


,n-l 


each cont.aiiihig 9 control points (Soo Figure 5.21) i.e., 


Bm • — [Bin, Of B„i^] , • ■ . ) ^^7)1, 71- 


(5.20) 


where 


B ■ •= \B^ ■ B'^ ■ B^ 


(5.21) 


and 


B'’ 


•— ^rnyjj ^ 


(5.22) 


The superscript T used in (5.20) and (5.21) denotes block transpose. The matrices 
Brn^x and Bm are related to each other by the following equation, 

Bm-ri = -^m+i Bm (5-23) 

where Mm+i ks a K x K matrix termed as (m+ l)th level subdivision matrix. Consid- 
ering the block form of Bm, i.e. Bm = [Bm,o, • • • , Bm,n-i?\ the matrix Mm+i can also 



be \vrit4,eM as a block matrix 


Al 


rn+1 


f (^Aii 1 1 ) 1,0 


n,0 


\ I ) 1 ,71- I \ 


where (iV/,„+i),j denotes tlu; (i, j)-th l)lo('k of and each block contains 

9x9 scalar elements. In particular we have 


■^rn+1,0 1 ) l, 0 '^m ,0 + (-^^171+1)1,11— 1 .^ 7 n,n—l 

-f? 77 i+l,l = (^^771+1)2,0^771,0 H + {Afm+l)2,n-lBTn,n-l 

~ (-^^771+1)1 , 71 - 1 *^ 771,0 d" (■^771+1)1,0-^771,1 d' ■ ■ ■ d" { A [^^ i ) i ^ ji — 2 Bm , n—l 

which follows from the nature of the subdivision scheme and the indexing scheme 
adopted for the control points. Therefore we liave 

(•^''■^ 777 + 1 ) 2,1 “ {Aim+l)l,0 

(^ 77 + 1 ) 2, 2 = {AI,n+\)\,\ 


(■^■^ 771 + 1 ) 2 , 71—1 (-^^771-) I ) l,ri-2 

(M77I+ 1)2,0 =" (-^^ 777 + 1)1, 77 - 1 - 

One can sec similar relation holding for consecutive rows of blocks. Therefore denoting 
(M,n+i)i,j by Mrl+i, j = 0,l,...,n~ 1, Alm+i acquires the block circulant form [25] 




m;’+, 

*Ci 

■■■ M,";b 

MZl\ 

Mi+i 

■■• Af";b 

VK+i 

Ml+i 

■■• 

by bcirc(M^+i, . . 



matrix A4,n+i it is enough to determine the block matrices A^m+i^ • " 1 ^m+i Note 
that these block matrices satisfy the equation 

•St77+1,0 = -^TTi+l-®”*!! d- 




(5,24) 



It is ('asy io observe (sec Figure 5.21) that 

Bmi lft = n + 6'Zll „_, + till „_1 


■'^''jn+1,0 ^ B„i^o "b q 4- sZJjij Q + Q 

B'm + lfi = ^ -^m.O + ^^m,0 + ^ 

-^m+l.O = ^'■^m.O + ^^m,0 + 

•^m+1,0 ~ ^^m,0 "b "b "b ^-^m.O 

C+1,0 = + sB^,„ + + tBi.,0 

Bf„4.t n = + sB,^„n + nB'n + n 


Vn,n“-1 


n5 

-^711+1,0 


n,0 + ^B, 


'm,0 + 5-^1, 0 + *-®m,0 


■^TK+l.O = + sBto,o + tBl, 

i4+1.0 = r<0 + 6<o + -<,l+^^7^ 

where = a[n,in + 1), h' = p{n,m, + 1), = 

/IM.?;?, 4- ll and t = '7('1,7n. + 1). 


■ + 


Comparing th(' above; system of eciuations with (5.24) we get. 




Let, n'i)ms<!iit„s Uic onUy of Uic matrix M,n+i. It is easily soon 

that all the eutiies of are non-negative and the following relations are satisfied: 

j=i cos2(/i/2"*) 

K (5.25) 

S"’ 


Now the patch correspotuling to the j th unit of ilm is written in vector form as 

:= := bm{u, v,j) Brn, (u, v) € w{Tn), (5.26) 

where 6rn(w,n,j) is a row vector of tensor product trigonometric B-spline functions 
k = 1/2,... K: 

b 7 n{n,v,j) := [bl,/'a,v,j),...,bl^/u,v,j)]. 

The bfXsis hniction.s b'^n{u,v,j), k = 1/2,..., K are localized and only 9 of them are 
nonzero in any unit of ll,,,. Denoting hy • - Mi and repeating 

(5.23) we get, 

• • • Ml Bo =: Bo. ' (5.27) 

Therefore, convergence of the subdivision scheme depends upon the eigen values of the 
matrices and hence on its factors M^+i- We study this aspect in the following 

subsection. 


5.6.2 Eigenvalue Analysis 

Recall that M,n+i = bcirc{Ml^^, ■ ■ ■ , AC;!). We now use the Fourier matrix method 
[25] to diagonalize the matrix Mm+i- 1'*^^ us define w := exp{i2'n/n) and 

M^+i =; ^ = 0, 1, . . . , n - 1. (5.28) 

J=0 

Moreover, for any discrete set {poi • • • tPn-i} define by 

jf = 

j=0 


(5.29) 



Note Unit if k - 1.2, ...9 arc eigenvectors of ; = 0, 1, . . . , (n/2), then 

following Corollary 1.5 the vectors j = 0, . . . , [n/2], t = 1, 2, . . . , 9 


where 

Hk •- [Vj,k, %k, ■ • . , Vj^kf 
form a complete set of eigenvectors of M^+i- 


(5.30) 


Let us (lefiiH' := 5" + b^ w~^ + 4 - /;2 ^ - Q, 1, . . . , n - 1. 


Then for A: = 0, i, . . . , 71 , — 1 we have 
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(5.31) 


The eigen values of are 

1 


b\ 


8 cos'^ ( ) cos ( ) 4 cos'^ ( ^ ) 16 cos^ ( ) cos^ ( ) 
Therefore by Corollary 1.5 the numbers 

i>', j = 0,l,...,[n/2] 


, 0 , 0 , 0 , 0,0 


and 


4 cos-^i^y 8 cos2(5;^r)cos(^)’ 16 cos2(^) cos^^) 
are eigenvalues of Af^i+i- All other eigenvalues ol Mm+i ^'I'e 0. A simple calculation 
shows that 

= 5" + 2cos(— )(6’ /c = 0,l,...,n- 1. (5.32) 

71 

Therefore, the largest eigenvalue of Mm+i is := 5° = TosW-)^ 
eigenvalue 

7 , 2 + cos(27r/7>.) L, ^ ( 533 ) 

^ " 8cos2(/i/2"‘+0cos(/i/2''^) 4 cos2(/7/2^+i) 



is tlu' sccotid larg('sl ('if-('nvalu(' of ;U„., , ror all ptovidod (T < /,./2"'+‘ < tt/T. This 
range has l)oen found computationally. 

Let. M = lini„i. t i- Ihci limit exists because the sequence of entries in M^+i 

coav(!ig(!s. ISotf^ that M also has a complete set of eigenvectors {c,V 2 ,V 3 , . . . ,Vk}. 

I he ('igeiiv<rlucs 1, A^, A;), . . . , A;^- of A4 are the limit of the eigenvalues of In 

j)aiticuiai, the hugest (ugimvahie of M is 1 which is associated with the eigen vector 
e = ( 1 , 1 ,...,!). Let V 2 and be the eigenvectors associated with the second largest 
eigenvalue X 2 = A;! of M. 

Below we find A 2 , V 2 and for the cases n = 3,5 and n = 6 . The eigenvectors are 
obtained using Corollary 1.5. 

11=3: 

In this ctise A 2 = A.i = 7/lG. Tlie corresponding eigenvectors V 2 and V 3 of M (up to 
t.hree d(;cimal j)laces) ar(i 

V2 = (0.107, 0.0G5, 0.527, 0.494, -0.032, 0.497, 1.088, 1.069, 1.0, 

- 0.149, -0.463, -0.735, -0.317, -0.849, -1.189, -1.518, -0.995, -0.5, 
0.042, 0.397, 0.208, -0.177, 0.882, 0.692, 0.43, -0.074, -0.5 )^ 
and 

U 3 = (0.110, 0.497, 0.545, 0.081, 0.999, 1.086, 1.124, 0.532, 0.0, 

0.037, -0.191, 0.184, 0.387, -0.528, -0.112, 0.380, 0.66, 0.866, -0.147, 

- 0.305, -0.729, -0.468, -0.472, -0.973, -1.505, -1.192, - 0.866 

n=5: 

In this case A 2 = A 3 = 0.5385. The corresponding eigenvectors V 2 and ^3 of M (up to 
three decimal places) are 

V2 = (0.319, 0.458, 0.764, 0.683, 0.592, 0.872, 1.15, 1.067, l.p, 

- 0.055, -0.344, -0.131, 0.147, -0.584, -0.380, -0.199, 0.046, 0.309, 

- 0.353, -0.671, -0.846, -0.592, -0.953, -1.107, -1.272, -1.039, -0.809, 

- 0.163, -0.070, -0.391, -0.513, -0.005, -0.304, -0.588, -0.687, -0.809, 
0.252, 0.627, 0.604, 0.274, 0.95, 0.919, 0.909, 0.614, 0.309 ) . 



I';, (U.Ki'J. U.-Ml. U.HB?, O.IHJT, U.80G, U.G83, U.o82, 0.299, 0.0, 

0..3r,.l. (I.OIKK 0.8.10, 0.67, 0.812, 1.040, 1.273, 1.107, 0.951, 

^ 0.347, -()..304, -0.040, 0.205, 0.385, 0.588, 

- 0.318, -0.082, -0.762, -0.455, -0.999, -1.065, -1.147, -0.869, -0.588, 

- 0.254, -0.278, -0.607, -0.629, -0.314, -0.618, -0.913, -0.922, -0.951 


n=6: 

In tins Ccvse = 5/8. The eigenvectors V 2 and vs of M are 

V2 = {0.378, 0.554, 0.803, 0.722, 0.710, 0.925, 1.143, 1.06, 1.0, 

0.054, -0.128, 0.115, 0.311, -0.254, -0.033, 0.164, 0.318, 0.5, 

- 0.324, -0.681, -0.688, -0.411, -0.965, -0.958, -0.98, -0.742, -0.5, 

- 0.378, -0.554,’ -0.803, -0.722, -0.710, -0.925, -1.143, -1.06, -1.0, 

- 0.054, 0.128, -0.115, -0.311, 0.254, 0.033, -0.164, -0.318, -0.5 
0.324, 0.681, 0.688, 0.411, 0.965, 0.95, 0.979, 0.742, 0.5 )‘ 


u;i = (0.156, 0.467, 0.331, 0.058, 0.704, 0.572, 0.470, 0.245, 0.0, 

0.405, 0.713, 0.861, 0.654, 0.967, 1.087, 1.225, 1.04, 0.866, 

0.25, 0.246, 0.530, 0.596, 0.263, 0.515, 0.755, 0.795, 0.866, 

- 0.156, -0.467, -0..330, -0.058, -0.704, -0.572, -0.470, -0.245, 0.0, 

- 0.405, -0.713, -0.861, -0.654, -0.967, -1.087, -1.225, -1.04, -0.866, 

- 0.25, -0.246, -0.530, -0.596, -0.263, -0.515, -0.755, -0.795, -0.866 )‘ 


5.7 Convergence Analysis 


Let u.s define tlie spjuie by 

:= {x := {xi, . ■ ■ (^^ii ^i) ^ !■ 

The elements of is considered as a matrix of order K x 3. 



l.H r, l.c ih,. roltniui vector in whose all coordinates are zero except the jth 
roordinate, It is well known that the set of vectors e^, j = l,2,...,K forms a basis for 
R^' • 5 a'! X ■" {./'i , . ■ . , f- 1 hen x can l)e exiirossed 


as 


X — Cl .i:j + , . . + Xf^', 

Define ilx||,, = inaxj. ii.r,j| vvlu're ||.r,lj = niax(l;r!|, |x?l, la;?!}. 


\\i' have ment ioned in the i)revious section that M has a complete set of eigenvectors 
e, r.’ 2 , . • . , vi(. Th(‘S(> eigcm vectors form a basis for Let us denote this basis by 

B := {ni, V2,...,vk} 

where Vi = c. I'lieu for x 6 there exist row vectors Pi € z = 1,2, 

such t.hat 

X = WiPi + . . . + ZJ/cP/v'. 

Define tin' injrm l|x||„ := maxKi<K \\pi\\ where |lpi|l = max{lp|l, jpfl, |p?l}. 

Now consider tire two sets of basis elements of {ef}f and {vi}f • Let T be 

the transformation matrix corresponding to this change of basis i.e, 

where f’tOi = EHi Clearly T is an invertible KxK matrix and 

llxll„ < llTll 11x11, and HxH, < 1|T-'|| HxH,. (5.34) 


A secnieiu'c r, , rj, • • • , in R^' is said l.o be order o(r) if 


lim 

m->c»o 


kmlL 

im 


= 0 . 


The above sequence is said to be order 0(/’") it 


lim 

m-400 



< OO. 


Let X = Bo- We show that tlic sequence of configurations Afl^'x converges to a lim- 
iting configuration. First we break the operator AfW into three parts, one stationary 

and two nonstatiouary. 



Ivt II, s ilriiiic .s;,, A/,„ A/, 

'){/), 7?/) > all the cnt-iii's of is 
in M- Thrn'fuii' all fhc entries in S,„ 


III (. binee o'(n,7n) > a, /j(n, in) > /i and 
gieatei than or equal to the corresponding entries 
are non-negative. The following recursive relation 


for holds. 


Lemma 5.1. Foj- ni > 2 

ni- 1 

j = l 


Proof: Wc jirove the lemma by induction on m. 

Note that il/dl = A/] = M + S]. Thus 

- A/2A/<‘> = (M + S2)(M + Si) = M^-hMSi + S2Mi. 

This provi's (5.30) for in = 2. 

A,s.sunu> that (5.35) holds for ni. Now 

A/h'*' » = A/„,h = (M 4- 5,„+i)A<f("‘) 

‘ in-1 

= A/{ A'/"‘ 4 } 4 S,n+i 

j = l 

m-l 

= A/”‘+' 4 ^ A'/”‘+'-^' Sj + MSm 4 Sm+i 

j=i 

Therefore tlie lemma holds for r/i 4 1. This proves the lemma. □ 


Let us denote y„i := 13jLV ^ 

To show the convergence of the sequence to a limiting configuration in 

j^/cxs show the convergence of the sequences {Af'”x}, {ym} 
has been shown by Reif [63] that 

M^x = epx 4 Ar;(u2P2 + vm) 4 o(A^) = epx 4 o(l) (5.36) 

where A 2 is the second largest eigenvalue and U 2 and V 3 are the corresponding eigen- 
vectors of M. Therefore {A-f'”x} converges. 



'Fo sinnv till- of {y„i} and we need the following esti- 

inatf’S, whirli we (Iciivi* in the following leininas. 


Lemma 5.2. For tniij x € !R ' ’ ’ ami iii > 1, Ihr fallowing holds: 


(5.37) 


Proof: Recall that 

Let n.s denote the (?, j) th element of the matrices MrX and Mr by (Mrx)ij and {Mr)i,j 
rcsj)ect ively. I hen, 

K K 

!(A/rX).J < \Y^{Mr)i,,xi\ <Y,\Wi,s\ Ix’il. 

(ir::! J)=l 

Since |.r'l < |lxl|^ we luive 

K 

\{MrX)J<Y^\{Mr)i,s\ llxlle- 

. 1=1 


JU cosHh/2'') w._To 

So |(M.x)jV^-«^ 11x11, i = 1,2,3; i = 1, 2, . . . Therefore, 


K. 


P^rXlle < cos2(/i/2'-1) 


cos2(h/2’-)_|, ,, 


( 5 . 38 ) 


lienee (. 5 . 37 ) follows. O 


Taking limit in (5.38) 
Lemma 5.3. For any x 


as r -4 oo we get the following result, 
g jfg^xxs fii>l we have 


ljAtf"*x|i, < IML- 


(5.39) 


Now wc obtain an cstinnatc for ||S„.xl|, in tbo following lo.nn.a. 



Leniina 5.4. t'or ni > 1 mid mui x G wr. have. ■ 

llSmxll, < ^ ||x||, (5.40) 

for smiii l■on^iilnl ( '{li) luih'pnideiit of in. Hc.nce the sequence {Smx} converges to 0. 


Proof: L<*i denote ihe (?',j)th entry of Sm- Then for i = 1,2, and 

j ~ 1,2, 3 wo have 

K K 

j(.s„x),.,j i 52 i((.Au,„, - (j/),,,) .Ti;i < EiwJm - {Mk.\ 


jfrr 1 


X 




!< K 


.r-=l 


5=1 

^K 


Since 


> (-WloA- fot- ail i,k = 1,. Jv. Further E^=i(^)m = 1 ^nd 

<>oh2(/(./2"0 


K 


5=1 


Tlierefort% 


l(S’„,x),.jl < 11x11, (cos 2 ( 5 ^ “ V - 
where C(/i) = □ 


X ^ 3/1^ C{h) 


4"* cos2(/i) 4"^ 


In t.ho following leniina wc show that the the sequence (ym) converges. 
Lemma 5.5. There exists unique q € such that 

lim Ym = eg. 

771-400 

Besides, l|y„, - cg|l^ < ^ for some constant Ci- 


(5.41) 


Proof: First of all we show that the sequence {ym} is Cauchy. 



1-ui 711, li 1 w<» liavi» 


ffi * k ■ I 


' ■ _ m~i 

y,:..: V„. L '.S'jA;" V«™ Js Ji/U-I)x 

J-! ,-=l 

?ri I 

H (a/"’+*=- ■’ Sj Af'J-i) X - Sj x) 

?» 4 4' i 


;--rri 


TiK'refort’ 


l>l-- 1 

||yn, » * - y«. II, < ^ Sj X - Af"‘-^' Sj X 

in. 4 A’-* 1 

+ E ||'W""'‘-'S,M«-‘>x||, 


(5.42) 


For any ilnarc of j, donoit’ v := Sj M^^~^'^x. So v = ui /9i + D 2 A '^k Pk for 

soinr /f, L E\ ? =* 1, ■ ■ • , K. I’hea l)y Lc’mma 5.2 and Lemma 5.4 we have 


MUsi&INL- 

Moreovei" 

Jipn+k~] y _ yj 0^ ^ ^2 (02 H 1" t'K Pk 

and A/”‘“-' v = i-, /?i + A!]*"-' U 2 + • ■ • + A'^"^' Vk Pk- Therefore 
||A/”*+'=“-'v- A/"‘“-'v||^ 

= ||(A"‘ ' - Ar') V 2 P 2 + ... + PkI 

li A"‘ •' (A*, ~ 1 ) 1)2 p2 -I- . . . + A’r^' (Ax - i) Pk IL 


(5.43) 


(5.44) 


< A"'--'^ max lAj - 1| (JC- 1) ||r-'|| |1 t|| ||v||^. 
Since 0 < A^ < 1 for 2 < r < K and K is fixed we get 

y _ ^/rn-j y||^ < Cj l|v||^. 

For some constant C 2 . Thus by (5.43), we have 

11 m-7 II 

V - M™-' v|l^ < A, jrjjjjft) '1=^ 


(5.45) 



Furtli*'!. iisiiir, F''inni;i;, :>.2. aiul 5.4 wc get 


.’I/""''' •' ,s’ y*-' Ox < 

c. nr\c 


C{h) X 


coH^ih) Ai 


Suii.sliluliiig aiiti (5. 40) in (5. .12) we get 


* k >%?! S ^’{ 1 ^) X 


ni- 1 


J-i j=m 




Since 4^2 > 2 and the .sums 


j=j ¥ exist we get 

1 


ym-ik ~ ymllj, < C{h) ||x||^ (A 2 " + — ) 


(5.4G) 


(5.47) 


for ii generic- constani ( '(/() imh'peiulent of rn. This proves that the sequence {ym}TO =2 
is ('aiieiiy ami hence cunveigenl. 


Now, we show that (r).-ll) iioltis. Let us assume that \hnh~^ooym = y and y is 
expre.s.sed in fenns «»r eigenvectors of M as 

y = Vi (3[ + ?)2 -f . . . + Vk P'j^. (5.48) 

for .some /f,' L E*, t I, . . . , J\ . 'Fhen wc; have 

||A/ y - y|t < 11-'^/ y - A/ y„.||^ -t- p-/ y,„ - y„,.M||, + ||ym+i - y|L 
< i|y - ym\l + II A'/ s,n x||^ + ||ym+i - y||,. 

Therefore by Leinnnus 5.2 and 5.4 wc get 

II A/ y - y P < ||y - yml|, + 4 ^^^ II^IL + l|y-+i " ^ 11 - 

By taking in 00 wc* have' My = y. Therefore y is an eigenvector of M for the 

eigenvahu' I. 'I'hus y = vi ft[ := e(i[. Writing /ij = q, we get the required result. 

Since A 2 < 4 , t aking A: -4 00 in (5.47) wo get ||?/ - 2 /m|lg < Hence second part 
of the knnina follows. □ 


By (5.36), (5.-11) and the. fact that the sequences {5mM("*-^)x} converges to 0 we 
get 

Brn = e{Pi + ?) + A2'(^2P2 + vzPs) + o(A^) + 0(l)/4"*. (5.50) 

The main tluHU’em concerning the convergence of the subdivision scheme follows as a 
conseciuonce of the ahov(^ considerations. 



V'ltrlt ]>i liViI ij HU lit JlUiil /)!/ (fi.rtO). 


Then 


Proof: Let x,„ c- Thcit l>y (5.50) wc have 

J) ~ + f/) + o(l)|. (5.51) 

Let 7:,(.r; /i) l.c lh<* trigunoinelric B-spIino of degree 2 with knots {0, h, 2/i, 3/i} and let 
1\ dfiiote till' tianslatiiin 


/ ,(./'; li) ~ h){x *■ ?/t; /).), i 6 Z. 
It out he easily clH't kcd {see j5'ij) that 


i€Z 


os{h) ' 


(5.52) 


lit view of the local of hf„(», a,./) the sum i) ^ non-zero 

terms. 'I'he.se terms aie siipportcfi on J th unit of Si,,,. Furher, in view of the above 

equation 


hm{u,v,j)c 


T— s 


k-\ 


< 


cos^(/i/2’^) cos2(/i) 


(5.53) 


Ih'iice 5„,(i/. (/q i t/) I f)(l)| - nw^pl/iFq (/h +y)+‘^(f)- Thcrelore, x,„(u, u, j) -> 

(pi I q). 'This <'oitiplete,s the proof of the convergence. □ 


'Ihe following important con.Hequenee is immediate. 

P = PU (pi + <?) = U U (Pi + 9) 

m€Z 

i.e. adding the limiting value to tlic surface fills the hole of the surface and hence the 
limit stirface becomes continoiis. 

Below, we st udy the smoothness of this surface and show that the surface P is 
tangent plane cont.inous at p := pi + <1- 



5.8 .1 angciit l^laiie Continuity 


I„ rase nf statiunaiy sulHiivision sdioinos coniimiity of the tangent plane of the limit 
sui fares is slnwii 'fl.'lj with t lu' lu'lp of the characteristic maps associated with the 
Stilxiivisi.m srheines. Since, tJie subdivision scheme described above is nonstationary, 
the notion ol rhaiarteiistic maj) is modified to suit the present scheme. 


Definition 5.2. l’(»r the subdivision matrix M the sequences of maps 
'Ib,, ; Q,n ut > i thdined by 


'hm : ( W, e, j) -t V := bm{u,V,j) [V2,Vz], 

is calli’d a ,se(juence of characteristic mai)S, where V 2 and Vz are eigenvectors of M 
corresiionding to (!h' eigenvalue = A,-) and V is a /f x 2 matrix with the vectors V 2 
and a.s its cf^lumns. 


I’he rows of \ ' can be viewt’d as points in and .) can be considered as a 

two dimensional tensor product i.rigonometric spline surface with rows of V as control 
iroinfs. Note flint the vectins e -2 and Uz for the particular cases /r = 3, n = 5 and n = 6 
are obtained in Section 5.6. Tlie characteristic map for these cases are shown in 
Figure.s 5.22-5.2'}. Since at all the points on the map the vectors and have 
different, directions they are linear independent. Hence := det ^ ¥" 0 

for all (u. j) e S2„i- 

Thooreiii 5.T. The limit surface P has a continous tangent at p. 

Proof: By (5.50) W(‘ have 

Ihn - (/h t q) -t A^" {V 2 P 2 + + 0 (^ 2 ") + {o(Ar) + ^}- 

Therefore 

Xm(w, W,i) = 1;, j)|e (Pl + q) + ^2” ('^2P2 + ^iPi) + o('^2*) + ) + 4m 

= hmia, w, j)e (pi -f- q) + bm{% i) (’■’ 2 P 2 + Pz) 

+ v,i){o(A.^’‘) + 



Hut wi' i'V ! j .'kH r.j) (' 

r. j) V 


(i u 




0 ami 


'I'horel'oro 

Oh„i{u, v,j) e 
dv 


0 . 


Hciict' 

C I 

ihi 


x.uHi. r.j) 


y2 + ?’ 3 P 3 ) + ■^kn{u,V,j)^o{X^) + {o(A^) + ^^}}, 


0 ( 1 ). 


■% 

(Jv 


~ '^ 2 ” »'..;)(t-’ 2 P 2 + W 3 P 3 ) + j)|o(A^) + {o(A^) + 

and the <iim’tii)n of the normal vector to the surface x,„ for any value of {u,v) is given 
by the vecloi pruduct 


0 

Ou 


<1 

Ox' 


Xl"' I *’0')("2l^2 + vm) X ■^bm{u,V,j){V2P2+VzP3) 

0 nf ./vm\ , ^( 1)1 , ^( 1 ) 

•t- rr-lhn(xt,x! 

(/U 


0 


^ ^Jbnin,vJ){o(XT) + -^})- 


Ci,u.« and arc' boumh'd functions on any unit of Qrn and 4A2 > 1 for 

U'lliv*y. II II II 

any j ~ 1,2,..., A* we have 

4-Xm X Tf x,„{n, n,./) = x,„.„ x x„i,„{», u, j) 

frti 'fl'tf 

= Al'"|4"An(u, tUj)V2^J>m{u.vJ>2{P2 >< Pz) + ^bm{u,V, j)v2 — bm{u,V, j)v^{p2 X Pa) 


\ du 


0 

+ j)e 3 - 7 r Aii(u, n, 

t/u c/t' 


P 2 ) + ^b,n{u,vj)vz-^b^{u,V,j)V 3 {pz X P 3 )} 


+ 


(<.(.vp) I ()(!>; 


|rii / 

f? 0 


■^}U^i,vJ)V 2 ^-brn{u,V,j)v,) 


.. ( <An \ 

+ A./‘ (o(l) -f Q-„ d- 

= A:f |A„,(n,e,./)(p.iXp3) + a(l)} 

Moreover A,„ is n, s.ul ,H,si,,ivc tor all (u,». j). Therefore the norrrrali.od normal 

vector is givc'ii by 

“ iF;;;;;;('Vha >< ■w('v>.j)ii 

A,„ 1 1 1 . 11 . _> ^ ’’ -Y, ;= n W 

f A„.{ h, n, i) (P2 X Pa) + 


(5.54) 



I Ilf luti iua! i o t hn surface at the point p. □ 



Figure 5.23: 'I’lu* veetDi I' (diainamls) and the characteristic map forn = 5 



Figun' 5.24: 'J'he vector I' (diaHiomis) and the characteristic map 'I/o for n - 6 

5.9 Applications and Examples 

In this section \v(^ giv<* some examples illustrating oui subdivision scheme. One 
the advantages of t he scheme is t.hat taking different values for h we limit surface 
different sizes. 'Fhus h is used }is a design parameter. 
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Response to examiner’s queries 


Suggestion no 1: ” + Theorem 1.4 : Block-circulant matrices are not necessarily diagoniz- 
able In particular, the blocks themselves are not necessarily circulant. Hence Fj^AjP^ does 
not make sense. The best one can hope for is a transformation to block diagonal form, but 
this is sufficient for all application in scope.” 

Response: 

Theorem 1.4 (Diagonalization of a block circulant matrix) A block circulant matrix 
A — bcirc{Ao, Ai, , An-i) has the block diagonalization form 

A = (F„ X FmYdiagiFrr^AoF:;^, . . . x FJ. 

Suggestion no 2: ”+ Page 25 : This is my serious observation. The weights for computing 
the new points do not sum up to 1. Thus the rule is not affine invariant and therefore not 
meaningful in a goemetrical sense (results depend on the specified coordinate system). I am 
sure that the artifacts which are visible in Figures 5.25 through 5.29 are due to this problem. 
The choice of weights has to be corrected, and it should be explained how this choice is mo- 
tivated (something quite vague” natural generalization of the weights for the regular case” 
will be sufficient).” 


Response: 

This is a non-stationary scheme. Thus the sum of the mask computing newpoints OjS are no 
equal to one. Instead sum of the masks at mth level of iteration is equal to r^, = cos^{h/ 2 ^y 
The choice of weights for the masks is a natural generalization of the weights for the regule 
case. Since roVi-'-roo — we normalize our scheme by multiplying cos^(/i) to tl 

limit surfaces of the non-normalize scheme. As a result our scheme has now convex hr 


property and translation invariant property. Thus they are now useful in geometrical sens 
We are not normalizing the weights of the mask since some of the techniques of conv^ 
gence analysis, adopted in Section 5.7 may not follow for the normalized weights. We a 
describing the properties of limit surfaces of the normal scheme in detail in a separate sectic 
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1 Propert.ic's of limit surfaces 



Figure 1: Or^ml .vnirol mesh and three iterations of cur nonstationary suUivrswn scheme 
(normalized) for h -- O.tHIB. 



Figure 2: Convex Hull property, h - 0.698. 


, , ■ invariant To show this consider the initial control mesh ^ 
Our scheme is translation invar • , , i r written 

a-. , ,H ns discussed in Section 4. The initial set of control points Bo. wnt 

and Its regular part U . . ^ ^ 

in the form of a i.ioch rnatn., ^s. ^ 

control mesh by a veTtor t € R • . . « B 4. et corresponding to the new set 

[/; satisfies [/; = (4 + ° 

of control points. ^ of the subdivision scheme {Se} for the control points 

Let X and X' be the limit surfaces of he s ^ x+t. LetXfn e 

[7° and U^' respectively. It is observed m Sectio respectively 

,-m .. . r-m s. 

corresponding to a parameter value {u,v,j) 
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Since e is an eigenvector for all matrices Mm associated with the eigenvalues we 

get By (4.7) we have bm{u,v,j)e = for all {u,v,j) G 

Therefore 

, 1 

— Xm + 2(U\^‘ 

COS^{h) 

Since the limit surface P is continuous, the extraordinary points are also translated by 
cos%j ^' ^ whole, all points in the limit surface P are translated by vector ^ 3 ^ 1 ^ t and 
hence the normalized limit surface cos‘^{h)P is translated by the vector t. 

We also note that the normalized limit surface cos^(h) P lies in the convex hull of Bq (see 
Figure 2 ). 

In case control meshes haveing nonpositive solid angle the normalized limit surfaces show 
saddle point behaviour near extraordinary points. The Doo-Sabin subdivision surfaces in 
this case look better than our normalized limit surfaces (see Figure 3). But if all the vertices 
of the control meshes have positive solid angle, our normalized limit surface behave nicely (see 
Figure 4). 

Finally we observe that our scheme is a dual scheme. 



Figure 3: Comparison, (b) Doo-Sabin (c) Our scheme, h — 0.5. 
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Figure 4; Comparison, (b) our scheme, = 0.5 (c) Doo- Sabin. 

Suggestion no 3: ”4- above theorem 5.7 ; The statements of the linear independce of the 
vector {d/du Xm} and {d(dv Xm} is not proven. I understand that such a proof would be 
very laborious ( and perhaps not woth do be carried out in detail), but a figure showing the 
parameter lines of the characteristric map might be used to support the claim.” 


Response: 

Linear independence of d/du Xm and d/dv Xm is required for showing Am(u,v,j) := 
^ 0. A figure showing parameter lines of the characteristic map might sup 
port this claim. But we are providing figures which shows part of the surfaces Aq for th 
cases n = 3, 5 and n = 6 respectively. The surface Ao consist of 3n patches. The figur 
below shows 3nth, first and second patch of Aq for the cases n = 3, 5 and n = 6 respective!’ 
They are captioned as Jn. Clearly, these patches never intersect the plane z = 0. Becau? 
of block symmetry of the matrices Mi other patches of Aq behave identically, i.e., they a 
nonzero. 
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